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In  this  thesis  we  study  two  different  models  for  PDE-based  image  processing. 
Both  model  the  removal  of  noise,  also  referred  to  as  smoothing,  from  digital  images 
while  retaining  essential  features,  such  as  edges,  and  both  take  the  restored  image, 
represented  as  a function  defined  on  a rectangle  D C R",  to  be  the  solution  to  a 
minimization  problem  over  BV  space. 

The  first  model  uses  an  adaptive  total  variation  (ATV)  functional  defined  on 
BV  space.  We  first  define  the  ATV  functional  for  functions  that  are  not  necessarily 
in  any  Sobolev  space.  This  space  is  the  a-BV  space,  where  cv  is  a chosen  function 
to  locally  control  the  amount  of  smoothing.  Then  we  derive  important  approxima- 
tion and  compactness  theorems  concerning  functions  in  a-BV.  Having  defined  our 
functional  and  proven  existence  and  uniqueness  of  a solution,  we  then  study  the  asso- 
ciated time  evolution  problem.  Here  we  define  a weak  solution  u{x,  t)  to  this  problem 


VI 


and  prove  its  existence,  uniqueness,  stability,  and  asymptotic  behavior  as  t — > oo. 
We  prove  that  u{x,t)  weakly  converges  in  to  the  solution  u^o  of  the  original 

stationary  problem.  In  addition,  we  demonstrate  some  numerical  results  of  the  time 
evolution  ATV  model  as  well  as  prove  the  existence  of  a solution  for  an  updated 
ATV  functional.  Also  discussed  is  an  updated  version,  where  the  parameter  function 
a depends  on  the  solution  u and  not  on  initial  noisy  image. 

The  second  model  uses  a functional  which  smoothes  the  image  where  its  gra- 
dient norm  is  below  a certain  threshold  e,  that  is  where  |Vu|  < e,  using  either  the 
Laplacian  or  a regularized  p-Laplacian  for  1 < p < 2,  and  retains  edges  where  its 
gradient  norm  is  above  the  threshold  {\Vu\  > e).  We  in  fact  prove  that  the  solution 
u is  smooth  where  |Vw|  < e. 
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CHAPTER  1 
INTRODUCTION 

1.1  Purpose  of  Investigation 

As  mentioned,  this  dissertation  is  the  study  of  two  methods  of  digital  image 
processing  from  a partial  differential  equation  (PDE)  approach.  However,  instead  of 
working  in  the  more  traditional  context  of  Sobolev  spaces,  we  work  in  the  space  of 
functions  of  bounded  variation,  or  BV  space.  As  will  be  explained  in  the  first  chapter, 
BV  space  is  a more  natural  space  for  images,  represented  by  functions  defined  on 
some  rectangle  D C R",  to  belong.  Both  problems  studied  here  relate  to  image 
denoising,  that  is,  removing  noise  from  corrupted  images  while  retaining  essential 
features  of  the  image.  And  both  models  involved  remove  noise  while  retaining  edges 
in  the  image  by  utilizing  a built  in  ’’selective  smoothing”  feature. 

1.2  A Brief  Overview 

In  the  second  chapter,  we  introduce  the  total  variation  and  adaptive  total 
variation  models  for  image  denoising  and  some  of  the  relevant  results.  We  then 
define  a-BV  space  as  a generalization  oi  BV  space,  given  an  appropriate  function 
a.  Then  we  define  a minimization  problem  over  BV  using  the  a-BV  semi-norm, 
the  solution  to  which  is  the  restored  image.  We  then  investigate  an  associated  time 
evolution  problem  and  prove  existence,  uniqueness,  and  stability  of  the  solution  as 
well  as  show  that  its  asymptotic  time  limit  is  the  restored  image.  A few  numerical 
results  are  presented  for  this  model  as  well  as  a brief  investigation  of  another  adaptive 
total  variation  model,  this  time  using  an  ’’updated”  version  of  the  selection  parameter 
function  a. 


1 


2 


For  the  third  chapter  we  study  another  selective  smoothing  model,  this  one 
based  on  a model  for  plasticity.  As  in  the  first  problem  we  take  the  restored  image  to 
be  the  solution  to  a minimization  problem.  The  main  focus  of  this  chapter  is  proving 
partial  regularity  of  solutions.  We  first  consider  the  simpler  case  involving  the  2- 
Laplacian,  then  focus  on  the  model  using  a regularized  p-Laplacian.  Some  results 
from  the  2-Laplacian  case  are  also  used  for  the  p-Laplacian  case. 


In  order  to  discuss  problems  in  PDE-based  image  processing,  we  must  first 
include  the  necessary  Banach  spaces  and  notation  related  to  the  study  of  PDF’s. 
These  include  Sobolev  spaces.  Holder  spaces,  and  of  particular  importance  to  image 
processing,  the  space  of  functions  of  bounded  variation,  or  BV  space.  It  is  assumed 


to  the  above  spaces,  as  well  their  use  in  PDF  theory.  For  convenience  we  include  a 
brief  summary  of  these  spaces.  We  start  with  weak  derivatives  and  Sobolev  spaces. 
Let  0 C R"  be  open. 

Definition  1.3.1  If  u,v  E and  a = (ai, a„),  then  v is  defined  to  be  the 

a^^-weak  derivative  of  u,  written  v = if 


for  all  test  functions  (f  E C^. 

From  this,  the  Sobolev  kF*’^(D)  spaces  for  integer  k > 0 and  real  p > 1 are 
defined  by 

Definition  1.3.2 

= {u  : Q ^ R|  V|a|  < k,  D^u  exists  in  the  weak  sense  and  D°‘u  E L^{H)} 


1.3  Basic  Function  Spaces  and  Notation 


that  the  reader  is  familiar  with  the  embedding  and  compactness  theorems  as  relates 
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The  Sobolev  norm  is 

l|^||iy'=’P(n) 


and 


for  1 < p < oo 


||u||H^*,oc(n)  = |D“u|£,oc(n)  forp  = oo. 

|a|<fc 

We  note  for  the  special  case  of  p = 2 we  write  = W^’^{VL). 

Next  we  have  the  Holder  spaces.  Let  k be  any  nonnegative  integer,  0 < 7 < 1, 
and  a = (ai, ...,  o;„)  as  before.  Define 


I i(o)  _ I I 

1^1^^  =max|u| 


“In  = and 

\a\=j 


\U 


(7)  _ 


|u(a:) 


= sup 

x,y^Ct  1^ 


Then  we  say  u G if  u G and 

k 


Y + Y < °° 

j=0  |q|=A: 


with  norm 


1^11 


c*’T(n) 


Y + Y 


(7) 


j=0 


The  most  important  Banach  space  we  will  use  is  the  space  of  bounded  varia- 
tion, or  BV  space.  This  is  essentially  an  extension  of  BV  functions,  however, 
now  may  include  jumps  and  their  gradients  are  interpreted  as  measures. 


Definition  1.3.3  A function  u G V-{Q)  is  in  BV{Q)  if 


=(n) 


< 1 > < 00 


Mnv(Q)  - |y  ^ div{(f)dx  : p G Co(D,R"),  \p\ic 
with  norm  |iu||Li(n)  + Mijv(n)- 

We  also  write  the  seminorm  [u]BV(n)  as  / |Vu|.  This  is  the  notation  we  will 

Jq. 

More  will  be  said  about  the  gradient  measure  in  chapter  3. 


use. 


CHAPTER  2 

SELECTIVE  SMOOTHING  USING  A DAPTIVE  TOTAL  VARIATION 

2.1  Introduction 

This  chapter  is  an  investigation  of  the  mathematical  theory  for  adaptive  total 
variation  (ATV)  regularization,  a powerful  technique  in  edge  preserving  and  noise 
removal,  which  has  been  effectively  applied  to  image  restoration.  Let  uq,  defined 
on  a rectangle  12  C R^,  be  an  observed  image  which  is  the  result  of  the  true  image 
Uoriginai  with  added  noise  77,  i.e.. 


Our  goal  is  to  try  to  approximate  Uoriginai  as  best  as  possible  from  the  observed  image 

Uq. 

In  recent  years  much  work  in  image  restoration  has  been  done  using  total 
variation  (TV),  see  for  example  Rudin,  Osher,  and  Eatemi  [22],  and  the  results  have 
been  promising.  The  restored  image  is  then  taken  to  be  the  solution  to 


for  the  constrained  problem.  Here  a > 0,  a > 0 are  chosen  positive  parameters,  and 


for  the  unconstrained  problem,  and 


min  TV{u), 


U 


subject  to 


11“  - “o||i2(n)  = 


TV {u)  is  defined  by 


TV  {u)  = sup 


U 


(2.1) 
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for  any  function  u G L^{Q).  We  note  that  (2.1)  is  precisely  the  BV{Q.)  seminorm 
of  u,  as  stated  in  the  introduction,  if  the  above  is  finite.  This  definition  for  the 
TV  functional  does  not  require  differentiability  or  even  continuity  of  u.  In  fact  one 
of  the  main  advantages  of  using  TV  functional  for  image  restoration  is  that  jump 
discontinuities,  that  is  possible  edges,  are  allowed.  The  example  below  demonstrates 
this.  In  fact  the  only  solution  to  the  problem  below  is  one  with  edges. 

Example  2.1.1  Let  Q,  — (0,2)  and  define  the  functional  F on  LP‘{Ll)  by 

F{u)  = TV{u)-\-  I {u  — Uo)‘^dx  with 

Jq 

fO  ^/0<a;<l 
uo{x)  = I 

[2  ifl<x<2. 

Later  we  will  prove  that  there  is  a unique  solution  u to  this  problem  with  TV{u)  < 00. 
We  claim  that  the  solution  to  the  problem 


min  Flu) 

u6L2(Q) 


is  not  in  W^’^(f]). 

Ln  fact,  let  f G C(f2). 

Case  1 : There  exists  an  x ^ (0, 1)  such  that  f{x)  < 0 but  f{x)  < 2 for  all 
xe[l,2) 

a.  If  f{x)  < 0 for  all  x G (0, 1)  then  let 


[ sup^e(o,i)  /(^)  y e (0, 1), 

u{y)  = 


I sup^e[i,2)/(a;)  ye  [1,2) 

if  supj.g(o,i)  fi^)  7^  sup^-gji  2)  /(^)-  Then  we  can  see  that  F{f)  > F{u).  If 


sup  f{x)  i=-  sup  fix) 
xe(o,i)  xe[i,2) 


set  u{y)  = 0 and  apply  Case  4. 


= z<^qz  + qf-f  = ^q  + ^v^  {n)j  < (/)^ 


usy;  q-’z,  = vfj  -(n)^  < {f)j  fijMsp  udyx  ■q-Z^'ofi 

■(s‘i]  q-z  ] 

\ = 

‘(t‘0)  3 o J 

puv 


‘{{^)f-z) 


uira 


‘(^)/ 


(i‘o) 

uim 


q 

V 


P7  -daoqv  a?/^  Suififsi^vs  (g‘o)  uo  snonuiiuoo  sq  / p/  og 
'(S‘rJ  uo  z>  {X)f  puv  (r‘o)  UO  0^  (x);  f ssbq 

■n  = / 0?  f 9SVJ  fijddv 

puv  2 = n 79/  ‘z  = //  '{n)j  < (/)t-  7»?/7  99s  9ffi  -g  ^ /? 

(g‘T]  3 ^ 2 1 

I = {fi)n 

‘(l‘0)3/i  (x)/ (^‘o)9^jui  j 

9UJ/9P  (t  ‘0)  3 X IP  jof  0 < (x);  f I -q 


•(™)^ 

< {f)j  f^ixvvp  dm  udijx  Pt  '(l  ‘o)  3 x diuos  jo/  q = (x)/  fj  ts 

■z  < {x)I 

7077  yons  (g‘x]  3 X uv  spixd  9^977  puv  (x  ‘o)  3 x gv  jo/  0 < (x)/  : £ asBQ 

< U)d  os  'On  = n gyv;  uoy^ 


'Z  < {x)I 

7»?/7  yons  {z'\]  3 x puv  q > [x)/  py;  yons  (x ‘o)  9 x pixd  dJdyx  ■ Z 9sbq 

•n  = / o'l  f dsvQ  dfddv  puv 

0 = n 79s  0 = (x)/ //  '{n)j  < {/)j  davy  dm  uivBy  -g  ^ (x)/  /i 

(s‘l]3/i  (x)/ 'I 

} = 


‘(t‘0)  3/1  0 J 

797  udyjp  -Q  9 *x  dmos  jo/  q = (*x)/  udy^  '(x  ‘o)  9 x diuos  jo/  q < {x)/ /j  -q 
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Thus  we  conclude  that  if  f is  any  continuous  function  on  (0,2)  then  there  exists  a 
function  u with  a jump  discontinuity  such  that  F{f)  > F{u).  Thus  by  uniqueness, 
no  f ^ C{fl)  can  be  a minimizer  of  F. 


argument  f and  hence  f can  not  be  a minimizer  for  F. 

As  an  improvement  on  the  TV  functional,  Strong  and  Chan  [23,  9]  introduced 
the  weighted  TV  functional 


for  spatially  adaptive  (selective)  image  restoration.  The  function  a is  chosen  so  that 
a is  larger  away  from  possible  edges  and  smaller  near  a likely  edge.  Hence  we  allow  for 
greater  smoothing  away  from  edges  and  less  smoothing  at  the  edges.  Certain  choices 
of  a{x)  were  given  by  Strong  and  Chan  [23,  9],  and  their  numerical  results  were  very 
promising.  However,  many  theoretical  questions  such  as  the  existence  and  uniqueness 
for  the  minimizer  of  the  weighted  TV  norm  with  some  penalized  term,  and  for  the 
related  evolution  problems  when  a is  a function  on  VL  (not  only  piecewise  constant) 
remain.  The  question  whether  the  solution  of  the  evolution  equation  converges  to 
the  minimizer  or  not  as  t — > oo  also  remains  open.  The  goal  of  this  chapter  is  to 
investigate  these  problems.  Here  we  would  like  to  point  out  in  the  case  of  constant 
a,  results  were  obtained  by  Chambolle  and  Lions  [8]  for  the  minimization  problem 


Now  if  f E PT^’^(f2),  then  by  the  Sobolev  embedding  theorem  f is  a.e.  equal 
to  an  absolutely  continuous  function  f and  hence  F{f)  = F{f).  Thus  by  the  above 


minimize 


with 


where  cr^  is  known  and  A is  a continuous  linear  operator  on  IF{Tl),  and  by  Acar  and 
Vogel  [3]  for  the  problem 
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where  A is  a linear  operator  on  L^(fi)  and 


Jn(u)  = [ ^\Vu\‘  + D 
Jn 


is  defined  on  BV  space.  Results  were  also  obtained,  for  instance,  by  Vese  [29]  for  the 
functional 


and  its  corresponding  flow.  Here  a > 0 is  constant  and  (/?  : R — )■  R+  is  a convex, 
even  function  nondecreasing  in  R+  with  linear  growth,  and  K ; is 

a linear,  continuous,  injective  operator.  However,  the  results  for  the  flow  are  only 
in  the  dimensions  one  and  two,  due  to  the  methods  employed  there,  that  is,  using 
general  results  on  maximal  monotone  operators  and  evolution  operators  in  Hilbert 
spaces. 

In  addition,  Hardt  and  Zhou  [18]  consider  the  flow  related  to 


with  Dirichlet  boundary  data  u = g on  dVt  for  a bounded  domain  fl,u  = Uq  on  flx{0}, 
and  any  convex  linear  growth  functional  cp.  In  their  analysis,  they  approximated  the 
above  solution  by  the  flow  associated  with  = rje  * v{p)  where  is  the  usual 

mollifier  on  R”. 

In  this  work  we  shall  extend  the  results  of  Chambolle  and  Lions,  Acar  and 
Vogel,  as  well  as  those  of  Hardt  and  Zhou,  to  the  adaptive  TV  scheme.  In  particular, 
we  shall  develop  mathematical  theories  for  the  problem  of 


that  is 


= divxP>p{Vu),  {(fp  = V(/?) 


minf  («)  = TV^  + f 11“  - “ollii(n) 


(2.2) 
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over  BV{Vl),  and  its  corresponding  evolution  equation, 

dvi 

— = diw:,{a{x)ipp{Vu))  - P{u-uq) 
du 

— = 0 on  dO,  X Rt’, 
on 

u{0)  = uq  on  0, 


on  n X Up, 


(2.3) 

(2.4) 

(2.5) 


where  0 is  an  open,  bounded  Lipschitz  domain,  = [0, T],  ,5  > 0 is  a parameter, 
(p{p)  = \p\  on  R”,  (pp{p)  = Vp{p)  = p/|p|  on  R",  and  a{x)  is  a parameter  (smooth) 
function  used  for  edge  detection  and  to  control  the  speed  of  smoothing.  In  image 
processing  problems  a is  often  chosen  as 


a{x) 


1 

1 + k\WG(j  * Uq\^ 


(2.6) 


k being  a parameter,  and  being  the  Gaussian  filter  with  parameter  a.  The 
definition  of  a weak  solution  to  (2.3)-  (2.5)  is  similar  to  that  in  Zhou  [30]  or  Hardt 
and  Zhou  [18].  However,  our  method  of  proving  existence  differs  from  Hardt  and 
Zhou  in  two  ways,  the  first  being  the  approximation  equation.  Second,  and  most 
importantly,  is  the  way  most  of  the  necessary  estimates  of  the  approximate  solution 
are  obtained.  Finally,  our  use  of  an  approximate  PDF  to  the  corresponding  flow  is 
in  contrast  to  Vese.  Also,  our  result  holds  for  all  n. 


2.2  Preliminaries  and  Definitions 

We  start  with  the  definition  of  / a\Vf  \ for  functions  in  / G L^(fi). 

Jn 

Definition  2.2.1  Let  LI  be  a bounded  open  subset  of  R".  Let  f be  a real  valued 
function  on  such  that  f G L^{L1).  Also  let  a{x)  > 0 be  a continuous  real  valued 
function  on  O.  Then  we  define  the  a-total  variation  of  f or  a-TV  of  f to  be 

/ a\Vf\=  sup  < / f div{(/))dx  : \(f{x)\ < a{x)  \/x  e Ll\ 

Jn  0eci(n,R")  Un  J 

where  4>  is  a vector-valued  function  <p=((j)i,...,(j)ji)- 


10 


Definition  2.2.2  We  define  f G to  be  in  a-BV  if 


sup  W / div{(j))dx  : \4>{x)\  < a{x)  Vx  G Q > < oo. 

0eci(f2,R")  1 7a  J 

Definition  2.2.3  If  f & a-BV  we  define  the  a-BV  seminorm  by 

/ a|V/|  = sup  iff  div{(t))dx  : \4>{x)\  < a{x)  Vx  G 
7a  1 7a 


and  the  a-BV  norm  to  be 


a-BV 


= f o|V/| 
7a 


+ 


z-i(a)- 


In  the  sequel  we  will  write  the  above  norm  as 


Remark  2.2.4  It  is  easy  to  show  that  if  f E 

[ a\Vf\  = [ a\Vf\dx. 

7a  7a 

Remark  2.2.5  Note  that  if  f E BV{Q)  and  functions  a and  fi  both  satisfy  the  con- 
ditions of  Definition  2.2.1  where  a{x)  < fi{x)  for  every  x E fl  and  f E fi-BV,  then 
we  have  f E a-BV  and 

[ a|V/|<  f l3\Vf\. 

7a  7a 

This  follows  directly  from  the  Definition  2.2.3  since  |^/>(x)|  < a{x)  implies  |<?!»(x)|  < 
fi{x). 


The  theorem  below  establishes  lower  semicontinuity  for  the  o-BV  seminorm. 

Theorem  2.2.6  //{/,}  C and  f E L^(H)  is  such  that  fj  f in  Lf{Vt),  then 

/ o|V/|  < lim  inf  [ a\V fA. 

7a  7a 

Proof:  Let  (j)  E (70(11,  R”)  be  a vector  valued  function  such  that  |0(x)|  < a{x)  Vx  G 
11,  then 

/ / div(^)dx  = lim  / /,  div(^)dx  < lim  inf  [ olV/d. 

7a  J-^°°7a  7a 
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Now  take  the  supremum  over  (j)  to  get 

/ ci|V/|  < liminf  [ a\Vfj\.  □ 

Jn  Jn 

Next  is  an  important  approximation  result  for  functions  in  BV. 

Theorem  2.2.7  Let  f € BV  where  a{x)  > 0 is  continuous  on  fl.  Then  there  is  a 
sequence  {fj}  of  functions  from  C°°{Ll)  such  that 


lim 


[ \ fj  - f\dx  = 
Jn 


and 


lim  f o;|V/,|dx  = f q:|V/|. 

j^co 


Proof:  We  essentially  apply  the  argument  of  Giusti  [13]  with  an  important  modifi- 
cation. Given  e > 0 we  construct  the  covering  {Ai}  of  where 

Ai  A\  = O2, 


with 


Ojt  = G : dist{x,  dQ)  > = 0, 1,  2, 

and  where  m is  large  enough  such  that 


Next  we  construct  the  sequence  {/J  so  that 


(2.7) 


i=l 

where  rj^.  is  the  usual  mollifier  on  R"  and  {f>i}  is  a partition  of  unity  subordinate  to 
{Ai}.  We  then  choose  the  efs  such  that  the  following  four  conditions  hold  simulta- 
neously for  each  ef 


1.  Ci  < e 
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ix  < e2' 


2-  [ IVci  * - f^i\da 

JQ 

3-  [ \Ve.*{fV(^>^)-fV(|>,\dx<e2-^ 
Jii 

4.  support  r/,.  * {f(j)i)  C ^i+2  ~ ^i-2- 


Summing  over  all  i we  get 

p OO 

/ I/e  - /Ma:  <J2  ^ 

an  JO 


thus 


/e  ->  / in  L^(Q), 


and  by  Theorem  2.2.6, 

/ q;| V/l  < lim inf  [ o;|V/e|dx. 

Jn  ^-*0  Jn 

Now  let  g G Co(Q,R")  be  such  that  |^(x)|  < a{x)  Vx  G Then 

n OO  n 00  n 

/ /ediv(i;)da;  = V / (^e^  * (/0i))div(i/)dx  = V]  / /</>idiv(r/,.  * 
an  an  an 


(2.8) 


so 


/ /ediv(£/)dx  = / /div(0i?7,,  * + Y]  / /div((;^i77,;  * g)dx  (2.9) 

an  an  an 

OO  ^ 

- X]  / - f^(l>i))dx. 

Z=1 

Denote  the  three  terms  on  the  right  side  of  (2.9)  by  I,  II,  and  III  respectively.  Note 
III<  ellallioo  by  our  choice  of  the  e^s. 

By  uniform  continuity  of  a,  there  exists  an  increasing  function  ca  such  that 
uj{r)  — > 0 as  r 0 and  |o:(2)  — Oi{z')\  < cu(r)  for  all  z,  z'  such  that  \z  — z'\  < r. 
Consequently  a{z)  < u>{ei)  + a{x)  for  all  |a;  — < e,.  Now  write  g = ag'  where  g'  = Q 
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if  q;  = 0 and  \g'\  < 1 Then  for  i > 1 and  any  x e 

\(i>i{x){r]i*  g){x)\  = \(j)i{x){qi  * ag'){x)\ 

< / (j)i{x)gi{x  - z)a{z)\g' {z)\dz 
Ju 

< / 4>i{x)gi{x  - z){a{x)  uj{ti))\g' {z)\dz 
Jn 

= (j)i{x)a{x)  / gi{x  - z)\g' {z)\dz 
Jn 

+(j)i{x)uj{ei)  / gi{x  - z)\g' {z)\dz 
Jn 

< a{x)+uj{e). 

So  we  get,  for  z = 1, 

1=  f * g)dx  < f a\Vf\+u{ei)  / |V/| 

Jn  Jn  Jn 

and  also 

11  = ^ / fdiv{(j)ig,.  * g)dx  < 3 / a\Vf\  + 3oj{e)  [ \Vf\ 
i=2  ^ Jn—no  Jn 

< 3e  + 3u{e)  [ |V/|, 

Jn 

with  the  last  inequality  following  from  (2.7).  Therefore 

[ fcdW{g)<  [ a|V/|+a;(ei)  [ \V f\  + 3e  + 3oo{e)  [ |V/|  + 

Jn  Jn  Jn 


q\a\Loc 


for  every  g € with  |5'(x)|  < a{x)  Vx  G fl.  Hence  taking  the  supremum 

over  g and  then  taking  the  limsup  as  e ^ 0 we  get 


limsup  / a|V/(:|da;  < f a|V/|. 
0 Jn  Jn 

From  (2.8)  and  (2.10)  one  finally  has 

J|m  f a\Vf,\dx  = j a|V/|.  □ 


(2.10) 


Remark  2.2.8  If  f E BV{n)  n a G C(fi),  and  dQ  Lipschitz,  then  there 

exists  a sequence  of  functions  {/«}  C such  that 


(2.11) 


And  if  f E L°°{Q),  we  also  have 


ll/n||L°°(n)  < C'(^)||/|U°°(n)- 


(2.12) 


In  fact,  in  the  proof  of  Theorem  2.2.7  we  choose  the  e’s  to  satisfy  l-f  and  in  addition 


\ve,  * if  (pi)  - f(pi\^dx  < e2  \ Then  we  can  take  fj  E C°°(fl)  n VK^’^(Jl)  n L^(fi), 


Since  C°°(Sl)  is  dense  in  n L^(f2)  there  exists  for  each  fj  a sequence  {fj,k}  € 

C°°{Sl)  such  that  as  k ^ 0 


Then  (2.11)  follows  from  (2.12>),  (2.1A)  and  a standard  diagonal  argument  applied  to 
{fj,k}  to  obtain  the  desired  sequence  {fn}-  By  the  construction  of  {fj}  and  [fj^k]  if 
in  addition  f E we  obtain  (2.12). 

Theorem  2.2.7  now  allows  us  to  prove  a compactness  theorem. 

Theorem  2.2.9  Let  {fj}  be  a bounded  sequence  in  a-BV  where  a E C(Q)  and 
in  addition  a{x)  > 5 > 0 ^x  E Q.  Also  assume  that  C R”  is  such  that  09. 
is  Lipschitz.  Then  there  is  a subsequence  of  {fj},  also  denoted  by  {fj},  and  an 
f E IT{9)  such  that  fj  f strongly  in  Lp{9)  where  1 < p < and  weakly  in 


Jn 

such  that 


(2.13) 


\\fj,k  ~ fj\\L^{n)  ->  0 and  \\fj^k  ~ /j||wci(n)  ->  0. 


(2.14) 


L—^{9). 


Proof:  Since  0 < 5 < a[x)  \/x  E 9 and  by  remark  2.2.5  we  have 
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Therefore 

f |V*|  < C. 

Jn 

Thus  fj  is  bounded  in  BV  norm  and  the  theorem  follows  from  the  compactness  result 
in  Giusti  [13].  □ 

2.3  Minimization  Problem 
We  now  consider  the  minimization  problem 

min  [ alVuj  f (u  — Uo)^dx.  (2-15) 

BV{n)nLHa)  Jn  2 

Here  uq  is  the  initial  noisy  image,  /3  is  a positive  parameter,  and  is  a bounded  open 
subset  of  R"  with  Lipschitz  boundary.  In  the  sequel  we  will  always  assume  that  a 
is  a smooth  function  satisfying  the  conditions  of  Theorem  2.2.9.  In  practice  we  can 
take  a{x)  in  (2.6)  for  some  cr  > 0 and  A;  > 0.  Assuming  uq  G we  see 

1 

“ 1 + C'||uo|licx>(n) 

The  constant  5 in  Theorem  2.2.9  can  therefore  be  chosen  as 

1 

1 + C'll“o|lioo(n) 

To  conclude  this  section  we  verify  that  (2.15)  does  have  a unique  minimizer. 

Lemma  2.3.1  The  functional  in  problem  (2.\b)  has  a unique  solution  in  BV{Q.)  n 
L2(Q). 


Proof:  Clearly,  the  functional  is  convex,  coercive  in  BV (fl)  nL^(fi)  and  by  Theorem 
2.2.6  is  lower  semicontinuous.  So  by  standard  results,  (2.15)  has  a solution  in  BV{Q,)C\ 
The  uniqueness  follows  from  strict  convexity  of  the  functional  in  (2.15).  □ 

2.4  Flow  Related  to  the  Minimization  Problem 


To  motivate  the  definition  of  a weak  solution  to  (2.3)-(2.5)  where  (p{p)  = p, 


P 


<Pp{p)  = j^,  and  Rp  = [0,  T]  we  follow  Zhou  [30].  Assume  the  solution  u is  sufficiently 
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smooth  to  justify  the  following  calculations  and  that  (p  is  as  mentioned  above.  For 
arbitrary  v E L^([0,  T];  we  multiply  the  equation  in  (2.3)-(2.5)  by  u — u to 

get,  after  integrating  by  parts  and  using  the  formula  p(p)  — p(g)  > pp(g)  ■ (p  — g) 
(due  to  the  convexity  of  p), 


/ u(v  — u)dx  + / ap(Vv)dx  > / ap{Vu)dx—  / P{u  — uq){v  — u)dx.  (2.16) 
Jci  Jo.  Jn  Jn 

Then  integrate  with  respect  to  t to  get 


nu{v  — u)dxdt  + / ap{V 

! Jo  Jn 


v)dxdt  > 


napiyu)dxdt  (2-17) 

! 

— I I /3{u  — Uq){v  — u)dxdt. 

Jo  Jn 


On  the  other  hand  if  (2.16)  holds,  by  selecting  v = u + X4>  for  (f)  E C^(Q)  we  get  that 
u is  a solution  to  (2.3)  in  the  sense  of  distributions. 

We  are  thus  led  to  the  following  definition  of  a weak  solution  to  (2.3)-(2.5): 


Definition  2.4.1  A function  u E L^([0,T];  BV(fi))  is  called  a weak  solution  of  (2.3)- 
(2.5)  if  ii  = dfU  E L‘^{Cl  x Rr),  u(0)  = uq,  and  u satisfies  (2.17)  for  every  v E 
L\[0,T]-BV{Q)),  a.e.  sG  [0,T]. 


Before  we  continue  we  list  some  properties  of  the  smooth  approximating  func- 
tion of  p, 

p^{p)  = 

defined  on  R”,  which  will  be  used  in  the  subsequent  discussion: 

1.  p^^{p)  is  convex  in  p, 


2.  p],{p)  - p>0  Vp, 

3.  p^  ^ p uniformly  with  respect  to  p as  e — )■  0. 

In  fact  0 < |<p^(p)  — <p(p)|  < e ^p  E R".  Here  we  use  the  notation  pp  = Vp^.  The 
proofs  of  the  above  properties  are  all  straightforward. 
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Consider  the  following  approximation  problem  of  (2.3)-(2.5): 


du 

m 

eAu  + div(o;(pp(Vu))  — /3(u  — Uq)  on  x R- 

1—1 

oo 

du 

dn 

0 on  dTl  X Rp, 

(2.19) 

u(0)  = 

Uq  on  n. 

(2.20) 

where 

Uq  G C°°(0)  with  Uq  uq  in 

(2.21) 

||'f^olU°°  — C'(fI)||uo||L°°) 

(2.22) 

1 a(f{V ul)  < C {fl)  f a(f{'Vuo). 
Jn  Jn 

(2.23) 

The  existence  of  Uq  is 

from  Remark  (2.2.8)  if  uq  G BV{Tt)  n L°°(fl). 

The  idea  is  to 

prove  an  existence  result  for  the  above  quasilinear  uniformly  parabolic  PDE,  obtain 
bounds  for  the  solution  independent  of  e and  6,  and  pass  the  limit  as  e — > 0 and 
5^0.  Indeed,  the  lemma  below  provides  an  existence  and  uniqueness  result  for 
(2.18)-(2.20). 


Lemma  2.4.2  The  approximation  problem  (2.\%)-(2.2Q)  admits  a unique  weak  solu- 
tion u^'^  where  G L°°([0,  T];  i7^(Q)),  G L^([0,  T];  and 

f f {u^’^)‘^dxdt f ^\Vu^’^{t)\‘^dx-\-  [ ap{Vu^'^{t))dx 
Jo  Ju  Jn  2 Jn 

< / + o;(p(Vuo)d2;  + e 

Jn  2 

for  a.e  t E [0,  T], 

Proof:  By  using  the  Galerkin  method,  the  fact  that  pp  is  a monotone  operator  ([6] 
and  [20])  we  have  a weak  solution  to  (2.18)-(2.20)  such  that  G L^([0,  Tj;  L^(fl)), 
^e,s  ^ L°°([o,T];/7i(fi)),  and 

^ J (u'’‘)‘‘dxdt  + j + aip'(Vu'^’\t))  + ^(u‘* -ulf^dx 


I 

Jn 


-|VuqP  + o(p^(VMo)dx. 
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Since  (p{p)  < ip^{p)  < (p{p)  + e for  any  p we  arrive  at 

f f {u^’^)'^dxdt+  f + f ap{Vu^'^{t))dx 

Jo  Jn  Jn  2 Jq 

— / 

Jq  2 

We  also  have  the  following  L°°(n)  bound  for  the  solution  to  (2.18)-(2.20) 
obtained  above: 

Theorem  2.4.3  Suppose  Uq  6 L°°{Q.)f]BV{fl)  andu^’^  is  a solution  of  initial  bound- 
ary value  problem  (2.lS)-(2.20).  Then  we  have 

||^^’‘^||L°°(nxRT)  ^ C'(^)||^o||L°°(n)- 

Proof:  Let  M = ||uolU°°(n)-  For  any  A > 0,  multiply  (2.18)  by  — M)+, 

where 


g At^£,(5  __  g At^£,5  _ > Q 

0 otherwise. 


and  integrate  over  fl  to  get 
du^’^ 


r du^>^  r 

J —e-^\e-^*u^’^-M)+dx  + eJ  Vu^’^  e~^^  ■ e~^^Vu^’^dx 

+ f app{Vu^'^)e~^*  ■ e~^^Vu^’^dx  + P f {u^’^  - uo)e~^\e~^^u"’^  - M)pdx  = 0. 
Jq  Jq 


Then  since  the  last  three  integrals  are  non-negative  we  see 

du^’^ 


Let 


Then 


f ^ — e — M)+dx  < 0. 

Jq 

nt)  = ljHe-V-MM^dx. 


e,6 


2 Jo  dt 


Fin 

Q dt 


)dx 


I 

JQ 


-At 


£.5 


du 


dt 


e ;r— (e  - M)pdx f -Xe  - M)+dx  < Q. 

Jq 
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Therefore  I{t)  > 0 is  decreasing  in  t with  /(O)  = 0.  Hence 

[ |(e-^V’-^-M)+|"dx  = 0 
Jn 

Vt,  and  then  C — a.e  on  Q,  VA  > 0 and  Vt  > 0.  Letting  A ^ 0 we 

obtain 

u^’\t)  < M = ||4l|L~(n)- 

Similarly,  > —M  = — ||uo||£<x.(n)  by  multiplying  (2.18)  by  — 

and  using  M = ||'ug||i^(n).  Thus  ||u^’‘^iU-(nxRT)  < IkolU-(n)  < C(fi)||'txo||L=c(n).  □ 
Before  we  state  the  existence  theorem  we  need  the  following  lemma: 

Lemma  2.4.4  Assume  the  weak  solutions  of  (2.1d>)-(2.2D)  have  uniformly 

bounded  L°°{Q  x Rj^)  and  L^{Q  x R^)  norms  in  e for  and  respectively. 
Then  there  is  a subsequence  of  still  denoted  by  such  that  as  e — )•  0 

1.  ^ h weakly  in  L^(f2  x R7-)  for  some  h 

2.  weakly  in  L‘^{Q,  x Rj’)  for  some 
where  h = ii^  and  u'^(O)  = Ug. 

For  the  proof  see  Zhou  [30]  or  Temam  [25]. 

Theorem  2.4.5  Suppose  uq  £ L°°{Q)  n BV{Q.),  then  there  exists  a unique  u € 
L^{[0,T]]  BV{fl)  n L°°(f2)),  ii  G x Rt),  and  u{0)  = uq  such  that  u satisfies 
(2.17)  for  a.e.  s G [0,T]  and  every  v G L^([0,T];  RH(fl)). 

Proof:  Let  be  the  solution  to  (2.18)-(2.20).  By  Lemma  2.4.2  and  Theorem  2.4.3 
it  satisfies 

||'t^^’‘^||L°°(nxR7’)  ^ ||^olU°°(n)  < C'(f2)||Mo||L°°(n)  (2.24) 
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and 


I \ii^’^\^dxdt  -f 

/ 

a^piy  u^'^)dx 

J UxRt’ 

Ju 

< 

1 

[ \Vul\‘^dx  + C 

'n 

< 

1 

[ \v4\'^dx  + c 

'u 

(2.25) 


/ a(p{Vul)d. 

Ju 


X + e 
X + e. 


By  (2.25) 


/ 


\u 


nxRx 


+ ^{Vu^'^))dxdt  < ||u^’'||By(nxR^)  < C|h-' 


with  C — C(q;,  Q,T).  So  for  fixed  5 > 0,  there  exists  a subsequence  of  such 

that  as  e ^ 0, 


> vr  strongly  in  L^(fi  x R^)  and  a.e  in  f]  x R-b  and 


u 


^ 


if  weakly  in  L^[Q,  x R^). 


(2.26) 


Notice  that  by  letting  e ^ 0 in  (2.24)  with  fixed  5 we  have 

||'w‘^||L°°(nxRT)  ^ C'(^)lko|U°°(n)) 


(2.27) 


since  from  (2.26)  a.e.  in  fi  x R^.  By  (2.26)  we  can  also  extract  a subse- 
quence, still  denoted  by  such  that 


u 


^.<5  ^ 


u"  strongly  in  L^{Q)  for  a.e.  t E [0,  T]. 


Also  notice  that 


as  well  since 


— > u'^  in  L^(Q  X R^)  as  e — >•  0 


(2.28) 


(2.29) 


/ 


\u^'^  — 


HxRt’ 


dxdt  < C(n)||uo||L°°(u)  f \u^’^  — u^\  dxdt 

t/s7  X Rt 


by  (2.24)  and  (2.27). 
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As  is  also  a weak  solution  to  (2.18)-(2.20)  we  have  as  in  the  motivation  of 
Definition  2.4.1 

nu'^’^{v  — u^’^)dxdt  + ^ f f \Vv\^dx+  f f aip'^{Vv)dxdt 

^ Jo  Jn  Jo  Jn 


> 


- / /*  \Vu^’^\‘^dxdt  + f f aip^{Vu^'^)dxdt 

^ Jo  Jn  Jo  Jn 


— f f {u^’^  — Uq){v  — u^’^)dxdt 

Jo  Jn 

> j I aip^ (Vu^’^)dxdt  — ^ f [ ~ ul){v  — u^'^)dxdt 

Jo  Jn  Jo  Jn 

for  all  V € L^([0,T];  Now  let  e ->  0 in  the  above  inequality  to  arrive  at 

nii^{v  — u^)dxdt  + f f a(p{Vv) 

! Jo  Jn 

>liminf^y  J a(p^(Vu^’^)dxdt  ~ (3  J J — Uq){v  — u^’^)dxdt^ 

> f f aip{Vu^)dt  — P f [ {u^  — Uq){v  — u^)dxdt. 

Jo  Jn  Jo  Jn 

for  all  V G L^([0,  T]; and  hence  also  for  all  v G L^([0,  T];  5V(D))  by  Proposi- 
tion 2.2.8.  Here  we  used  (2.26),  (2.28),  (2.29),  Theorem  2.2.6,  the  fact  that  <p^(p)  > 
(p{p)  for  all  p,  and  uniform  convergence  of  to  p.  This  shows  that  is  a weak 


solution  of  (2.3)-(2.5)  with  initial  data 


Un. 


Additionally  from  (2.28)  and  Theorem  2.2.6  it  follows  that 

/ p{Vu^)  < liminf  f p{Vu^’^)dx  for  a.e.  t G [0,T]. 

J Q,  t/ 


(2.30) 


Thus  letting  e — )•  0 in  (2.25), 

lii^pdxdt  + 


0 X Rj* 


f </?(Vu‘^)  <C  f <^(Vuo). 
Jn  Jn 


So  G L°°([0,  T];  BV{n)  n L~(Q))  and  e L\VL  x [0,  T]). 
Again  we  extract  a subsequence  in  5 to  get  as  5 — > 0 


u 


u 


u 


u in  L^(D  X Rr),  hence  in  LP{VL  x R7-)  from  (2.27), 
— >■  w in  L^{Q)  for  a.e.  t G [0,f],  and 
^ « in  L^(0  X Rr). 


(2.31) 
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Finally  pass  to  the  limit  as  0 in  the  inequality 

nii\v  - u^)dxdt  + f f a(p{Vv)  > f f aip{'Vu^)dt 
! Jo  Jn  Jo  Jn 


n{u^  — uo)(v  — u^)dxdt 

; 


to  get 


nii{v-u)dxdt+  f f a(p{Vv)  > f f a(p{'Vu)dt 
! Jo  Jn  Jo  Jn 

~P  I / {u  — Uq){v  — u)dxdt 

Jo  Jn 

for  all  V e L2([0,r];  W(f2)).  Therefore  we  have  the  existence  of  a weak  solu- 
tion u to  (2.3)-(2.5).  Using  (2.27),  (2.30),  and  (2.31)  we  see  as  before  that  u E 
L°°([0,T];BV(n)nL°°(Q)). 

To  prove  uniqueness,  consider  two  weak  solutions  Ui,  U2  to  (2.3)-(2.5)  with 
Ui(0)  = «2(0)  = Uq.  We  have  the  two  inequalities 

nu\{u2  - Ui)dxdt  + / / aip{Vu2)dt  > a(f{Vui)dt 

! Jo  Jn  Jo  Jn 

— / P{ui  — Uq){u2  — Ui)dxdt 

Jo  Jn 


and 


nU2{ui  - U2)dxdt  + f f a(p{Wui)dt  > f f aip{Vu2)dt 
'■  Jo  Jn  Jo  Jn 

- / /3{u2  ~ Uo){ui  - U2)dxdt. 

Jo  Jn 

Adding  the  above  inequalities  and  combining  we  get 

n{u2  — ui){ui  — U2)dxdt  > / / /3{ui  — U2)^dxdt. 

1 Jo  Jn 


And  hence 


giving 


[ ^ / (^1  - U2fdxdt  < - [ [ !J{ui  - U2f‘dxdt  < 0 
Jo  dt  Jq  Jq 


||ui(-,  s)  - U2{-,  s)||l2(0)  = 0 for  £-a.e.  s. 


Therefore  Ui  = U2-  D 


2.5  Stability  and  Asymptotic  Behavior 


Lemma  2.5.1  If  Ui  and  U2  are  two  weak  solutions  of  (2.3)  with  initial  data  uio  and 
U20,  then  for  C-a.e.  tE  [0,T], 

ll'^l  ~ U2\\L<^(n)  < ll^^io  — 'i^2o||L«>(n)- 
Proof:  Let  M = ||uio  — U2o||L°°(n)-  For  £-a.e.  t we  have 


/ Ui{v  — Ui)dx  ~\-  / aip{Vv)>  / aip{Vui)—  / ^{u\  — uiq){v  — Ui)dx  (2.32) 
Jri  Jn  Jn  Jn 


and 


/ U2{v  — U2)dx  + / acp{Vv)>  / aip{'Vu2)  — / I3{u  — U2o){v  — U2)dx.  (2.33) 
Jn  Jn  Jn  Jn 


Define 


= u{  — {u\  — ^2  — M)_|_ 

and 

= U2  + {u\  — U2  — M)+ 

where  u\  and  are  the  approximation  functions  from  Theorem  2.2.7.  Inserting 
and  w^  into  (2.32)  and  (2.33)  respectively  we  obtain 


/ ui{v^  — ui)dx+  / aifCVv^)  > / aip(Vui)  — / ^{ui  — uiq){v^  — ui)dx 
Jn  Jn  Jn  Jn 


and 


/ U2{w^  — U2)dx  + / ai^[Vw^)>  / a^p{Vu2)—  / ^{u2  — U2q){w^  — U2)dx. 
Jn  Jn  Jn  Jn 

Now  add  the  above  two  inequalities  to  get 

/ ui{v^  — ui)dx+  / U2{w^  — U2)dx  > / a^p{Vui)+  / a(p(Vu2) 

Jn  Jn  Jn  Jn 
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Observing  that 

/ aip{S/w^)  + / aip{Vv^)  = / a(f{Vu\)  + / o;(^(V'Ui) 

>/n  Jo.  Jn  Jn 

with  ^ V and  ^ to  in  L?{Q)  where 


V = U\  — {Ui  — U2  — M)^ 


and 


(2.34) 


(2.35) 


W = U2  + {Ui  - U2~  M)  + 

we  see  after  letting  e — >■  0 that 

I Ui{v  — Ui)dx  + / U2{w  — U2)dx  > — I P{ui  — Uio){v  — Ui)dx 

Jn  Jn  Jn 

- / P{U2  ~ U2q){w  - U2)dx. 

Jn 

But  the  right  hand  side  of  the  above  inequality  satisfies 

p / (ui  - tt2  + U20  - Uio){ui  -U2-  M)+dx  > 0. 

Jn 

Thus 

/ ui{v  — U\)dx  + / U2{w  — U2)dx  > 0. 

Jn  Jn 

Hence,  using  equalities  (2.34)  and  (2.35)  for  v and  to  in  the  above  and  combining, 

/ (ui  — U2){ui  — U2  — M)+dx  < 0 
Jn 

which  implies 

/ |(tii  — U2  — M)+pdx  < 0 

Jn 


A 

dt 


and  therefore 


/ |(tti  — ti2  — M)+|^dx  < / \{uio  — U20  — M)^\^dx  = 0 

Jn  Jn 

by  the  choice  of  M.  Similarly  U\  — U2>  —M.  □ 

To  conclude,  we  investigate  the  asymptotic  behavior  of  the  weak  solution  to 
(2.3)  by  showing  that  the  solution  converges  weakly  in  L^(fi)  and  strongly  in  L^(0) 
to  a minimizer  of  (2.15). 
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Remark  2.5.2  It  is  straightforward  to  show  that  if  inequality  (2.17)  is  satisfied  for 
all  V G L‘^{[0,T]\ BV{il)),  then  (2.1Q)  holds  for  all  v G L‘^{[0,T];  a.e.  t. 

Then  by  using  Young’s  inequality,  inequality  (2.1Q)  implies 


)dx 


[ f / (^ 

Jn  ^ Jn 


uo)‘^dx  > 


for  all  V G BV{fl),  a.e.  t. 


Finally  we  prove  the  following  theorem  concerning  the  asymptotic  convergence 
of  our  solution  u to  (2.3)  as  t — )■  oo: 

Theorem  2.5.3  The  solution  u to  (2.3)  weakly  converges  in  L?‘{VL)  to  a minimizer 
Uoo  of  (2.13). 

First,  let  F be  defined  on  BV{Q.)  n L^(n)  by 

F{u)  = f q;|Vu|  f{u  — Uo)^dx. 

Jn  2 Jq 

Second,  we  recall  the  definition  of  the  subdifferential  of  a proper  convex  functional 
G : H (—00,00]  where  H is  a Hilbert  space  with  inner  product  < •,  • >:  the 
subdifferential  of  G at  u,  written  dG[u],  is  defined  as  dG[u]  = [w  G H\G{v)  > 
G{u)+  < w,v  — u > Vu  G H}.  We  also  let  D{dG)  C H he  the  set  of  all  u with 
dG[u]  0.  Noting  that  the  above  F is  defined  on  a subspace  of  F^(0),  let  the  operator 
dF{u)  be  the  subdifferential  of  F at  u so  that  F{v)  > F{u)  + Jqw(v  — u)dx  Yw  G 
dF{u),  \/v  G BV{Q)  riL^(f2).  By  Brezis  [6],  dF  is  a maximal  monotone  operator  and 
by  the  above  remark  our  solution  u to  (2.3)  satisfies  for  a.e.  t 

Using  the  above,  we  can  prove  Theorem  2.5.3  by  the  following  lemma  ([7]): 

Lemma  2.5.4  Let  (/?://—)■  (—00,  +00]  be  a proper  lower  semicontinuous  convex 


function  which  assumes  a minimum  in  H . Then  for  any  Xq  G D{dp),  there  exists  a 
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unique  function  x : [0,  oo)  — > H which  is  absolutely  continuous  on  [5,  oo)  for  all  5 > 0 
and  which  satisfies 


x{t)  G D{d(p)  for  all  t > 0, 
x{t)  G —dcp{x{t))  a.e., 
x(0)  = Xq, 

and  iw-limt^oo  exists  and  is  a minimum  point  of  ip. 

To  prove  Theorem  2.5.3  we  take  H = p = F where, 

r F(u)  \iueBV{Q) 

F{u)  = { 

( oo  ue  L‘^{n)\BV{n) 

X — u,  and  Xq  = uq.  As  in  the  proof  of  Theorem  (2.2)  ([18]),  the  function  u : 
[0,oo)  ->  L^(fl)  is  absolutely  continuous  for  all  nonnegative  t.  So  by  Lemma  2.5.4 
we  have  u{t)  — ^ Uqo  weakly  in  L^(fl)  as  t — >■  oo  and  Uqo  a minimizer  of  F in  Z/^(f]). 
By  uniqueness,  Uoo  € 1/^(0)  fl  BV{Q).  Since  u{t)  is  uniformly  bounded  in  BV{^), 
we  may  also  conclude  that  any  sequence  {u(t„)}  has  a convergent  subsequence  still 
denoted  by  {u(t„)}  converging  to  u^o  strongly  in  L^(0).  Hence  u{t)  — >•  Uoo  strongly 
in  L^(fl).  The  theorem  is  now  proved. 

Note  in  fact  that  the  minimizer  u = Uqo  from  problem  (2.15)  is  actually  in 
To  see  this,  note  that  the  inequality  in  Lemma  2.4.2  and  the  proof  of  Theorem 
2.4.5  imply  that  ||u(t)||Bv(n)  < ^'(Q, (5)||uo||By{n)  and  ||w(t)||L°°(n)  < C'(f^)||^to||L°°{n) 
for  a.e  t G [0,Tj.  Hence  by  compactness,  we  can  extract  a subsequence  {-{/(tn)} 
such  that  u{tn)  — >■  u in  L^{fl)  and  u(t„)  — )■  u a.e.  on  fl.  Thus  u G L°°(fi)  with 
ll'a||L°°(n)  < C'(fl)||?^o|lL°°(n)- 
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2.6  Numerical  Results 


For  the  numerical  experiments  we  approximate  (2.3)  with 
du 


= div(o;(^p(Vu))  — P{u  — Ug)  in  x R^, 
= 0 on  do,  X Rt-, 
u(0)  = Uq  on  Q X {0}. 


dt 

du 

dn 


We  write  the  above  as 
du 


dt 


= adiv 


Vu 


\A  + |v 


+ Va 


Vti 


u\ 


Ve+|V 


p{u  - ul). 


Using  forward  time  differences  and  the  Neumann  boundary  condition  we  com- 
pute n = 1,2, . . . , N , N = number  of  iterations,  by 


= ul  + Ai6W,) 


where  uF  = uo{xi,yj),  and 


Q{u)  = a div 


Vu 


The  term 


y/e  + I Vu| 


Va  • 


Va  • 


Vu 


Vu 


i/e-b  |Vup 


- p{u-ul). 


V^e+TV 


u\ 


is  discretized  using  the  following  scheme  developed  by  Osher  and  Sethian  [21]  to 
permit  the  development  of  discontinuities  at  object  boundaries.  Here  we  let 


Aj.  Uij 


Ax  Uij 


Ay  UiJ 


AyUiJ 


2 ^ 

'^hj 

^ij  + 1 ~ 


2 
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Figure  2.1:  Original  256  x256  image 


Then  we  use 


(Va  • = max(Aj;Q;ij,  0) Uij  + min(Aa;Q;i,j,  0)A+Ujj 

+ max(Aj,Q;ij,  0)A~Uij  + min(Aj/0;jj,  0)A+Ujj 

as  in  Osher  and  Sethian  [21],  while  the  term  + [Vup  is  computed  using  the  central 
differences  AxUij  and  AyUij. 

For  the  term 

\^e+\Vu\\ 

we  use  the  scheme  from  Rudin,  Osher,  and  Fatemi  [22]  with  central  differences,  that 


IS 


div 


Vu 


y/e-h  |Vu|^ 


+a: 


These  figures  demonstrate  the  implementation  of  the  above  numerical  scheme 
using  256  by  256  gray  level  images,  with  pixel  values  0 to  255.  Figure  2.1  shows 
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hMiMmsiMms 


Pm%W¥ 


Figure  2.2:  Noisy  image  with  SNR=1:1 


Figure  2.3:  Restored  image  using  ATV  with  A:=0.001,  cr=0.5,  e=0.01,  /3=0.001,  300 
iterations 
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Figure  2.4:  Restored  image  using  TV  with  ^ — 0.001,  300  iterations 


Figure  2.5:  Noisy  image  at  j/=200 
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Figure  2.6:  Dashed  line  is  restored  image  from  Figure  1 using  ATV  at  y=200,  dotted 
line  is  plot  of  original  image  at  y=200 


Figure  2.7:  Dashed  line  is  restored  image  using  TV,  /3=0.001,  300  iterations  at  y=200, 
dotted  line  is  plot  of  original  image  at  y=200 
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Figure  2.8:  Original  256  x 256  image 


Figure  2.9:  Noisy  image,  SNR=1:1.5 
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Figure  2.10:  Restored  image  using  ATV  with  A;=0.005,  (7=0.7,  e=0.01,  ^=0.002,  300 
iterations 

the  original  image;  figure  2.2  shows  the  original  image  with  added  noise;  figure  2.3 
the  restored  image  using  the  above  ATV  scheme;  and  figure  2.4  the  restored  image 
using  the  above  scheme  with  a = 1 (restoration  using  TV).  The  next  three  figures 
examine  a cross  section  of  the  images  from  figures  2.1  through  2.4.  Figure  2.5  is  a 
cross  section  of  the  noisy  image  which  includes  the  bottom  tip  of  the  thin  ellipse. 
Figures  2.6  and  2.7  demonstrate  the  effectiveness  of  the  ATV  scheme  in  preserving 
the  tip  of  the  ellipse.  By  plotting  the  cross  sections  of  the  restored  images  and  the 
original  image,  we  see  that  the  ATV  scheme  is  better  able  to  preserve  the  tip  of  the 
ellipse.  Figures  2.8  through  2.10  also  demonstrate  the  effectiveness  of  the  method  to 
preserve  thin  lines  while  reducing  noise.  For  better  display,  we  threshold  the  restored 
images  to  be  between  the  gray  levels  0 and  255. 


2.7  Updated  Selective  Smoothing 

In  this  section  we  will  investigate  problem  (2.15)  with  our  choice  of  a being 

1 


a{x)  = 


l + k\V{G^*u)\^' 


(2.36) 
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So  in  the  time  dependent  case,  the  selective  smoothing  functional  now  uses  a con- 
tinuously updated  version  of  the  image  u.  Fix  a and  denote  by  G.  Assuming 
sufficient  smoothness  we  can  formally  derive  the  Euler-Lagrange  equation  for  (2.15) 
with  a as  above  to  be 


‘2^kGi,{y  - x)wi{y)\V  u\dy^  + div  (q;(^p(Vu))  - /3(u  - uq)  = 0 


with 


Wi  - 


Gi  * u 
1% 


2r^  n 

a Gi^*U,  ^ = 0 

an 


{l  + k\V{G*u)\^Y 

and  where  Gl^  is  the  partial  derivative  of  G with  respect  to  its  zth  argument.  The  first 
integral  on  the  left  hand  side  was  obtained  by  using  the  definition  of  convolution  and 
then  reversing  the  order  of  integration.  In  order  to  extend  the  solution  to  the  space 
BV{Q.)  n we  multiply  the  above  hy  v — u,  use  convexity  of  cp,  and  integrate  to 

get 

“^kGi^iy  - x)wi{y)\Vu\^  dx^  {v  - u)  + j ap{Vv) 

> / a(/?(Vu) -h  / P{u  — Uo){v  — u)dx. 

Jn  Jn 

Notice  that  in  the  first  integral  on  the  left  hand  side  the  coefficient  function  of  | Vu|  is 

not  necessarily  nonnegative  as  in  Definition  2.2.1.  However,  the  definition  of  / Of|Vu| 

Jo. 

for  arbitrary  a G C(D)  can  be  extended  to  be 


/ dlVul  = / d'^'lVul  — a |Vu| 
Jn  Jn  Jn 


where  dc^{x)  — max(d(a:),  0)  and  a (x)  = — min(o;(x),  0).  One  can  easily  verify  that 
if  u G then 

/ a\Vu\dx=  / oi^\Vu\dx—  / a~\Vu\dx. 

Jn  Jn  Jn 

Note  that  if  a is  not  assumed  to  a nonnegative  function  then  we  no  longer  have  lower 


semicontinuity  of 


[ d|V 

Jn 


u\.  From  this  definition  combined  with  Remark  2.2.5  we 


have  the  simple  result  below. 
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Theorem  2.7.1  If  a ^ C(fi)  n L°°(r2)  and  u G BV{Q),  then  / q;|Vm|  < oo  and 

Jn 


[ «|V 

Jn 


u\  < ||q;| 


L°°{n) 


f |V«|. 

Jn 


Instead  of  studying  the  complicated  flow  corresponding  to  the  above  Euler- 
Lagrange  equation,  we  may  instead  consider  a simplifled  version  of  this  flow  as  follows: 


du 

dt 

du 

dn 

u 


= diVi(o;(/Jp(Vu))  — P{u  — Uq)  on  x Rp 
= 0 on  50  X Rp 
= uo  on  0 X {0} 


with  a as  in  (2.36).  Note  that  a{x,  t)  now  depends  on  u{x,  t).  Although  the  deflnition 
of  a weak  solution  of  the  above  PDE  is  the  same  as  Deflnition  2.4.1,  the  dependence 
of  Of  on  t greatly  complicates  an  existence  argument.  So  here  we  will  only  prove 
the  existence  of  a minimizer  of  (2.15)  with  the  above  choice  of  a.  Uniqueness  is 
complicated  by  the  observation  that  F is  no  longer  a convex  function.  We  thus  only 
prove  an  existence  result  for  the  stationary  problem. 


Theorem  2.7.2  There  exists  a solution  belonging  to  BV (Tl)  C\  (il)  to  (2.1b)  where 
a is  given  by 

1 

l + k\V{G*u)\^' 

Proof:  Let 

F{u)  = I Q!(a:)|  V(u)|  + f {u  — uq)^  dx 
Jn  Jn 

and  let  {«„}  be  a minimizing  sequence  for  F.  Then  ||u„||p2(n)  < M,  M depending 
only  on  D.  Now  for  any  u G L^(Q)  we  have 


where 


a{x) 


1 

1 + |VG*u|2 


< 

< 
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with  C depending  only  on  0,  a. 
Denote 


1 

1 + |V(j  * Un\^ 

by  an{x).  Then  since  ||iin||z,2(n)  ^ M hy  the  above  we  have  o;„  > <5  > 0 for  some  5 
depending  only  on  0,  a.  Thus 


F{Un)>  / (5|Vu„l+  / {Un~  Uofdx  > 5 / iVu„ 
7n  JQ.  Jn 


Therefore  / |Vu„|  < M and  by  compactness  there  exists  a subsequence  which  we 

Jn 

still  denote  by  {u„}  and  a u in  BV{Q.)  fl  i/(D),  for  any  1 < p < 2,  such  that  Un  u 
in  U’{Q).  Let 

a{x) 

Since  — > u in  LJ{Q)  we  have 


1 

1 + 


/ G^.{x  - y)un{y)dy  ^ - y)u{y)dy 

Jn  Jn 

for  every  x E Q.  Hence  an{x)  o;(x)  and  |o;„(a:)|  < 1 for  every  x E G,.  So  o;„  — >•  o; 
in  I/^(D)  for  any  1 < p < oo  by  Lebesgue’s  Dominated  Convergence  Theorem.  Now 
computing 


dan 

dxi 


'iUo,,(x~v)uMdy  )(/.  Gx,xA^  - y)un{y)dy 
{l  + \VG*Un\^)^ 

^ Un  {x  - y)un{y)dy  )(/.  Gx2X\  {x  - y)un{y)dy 


2\2 


(1  + |VG*u„|^) 


gives 


dan 

dxi 


da 

dxi 


pointwise  on  Q,  and  similarily  for 


dan 
dx2  ' 


dan 

dxi 


< M' 


Since 
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for  i = 1,  2 and  M'  depending  only  on  G we  again  have  by  Lebesgue’s  Dominated 
Convergence  Theorem,  > a in  for  any  1 < g < oo. 

By  definition, 


/ anlVunl  = sup  < U„div(0)dx  : |0i(x)|  < an{x)  V(a;) 
Jn  <^>ec^(n)"  Un 


e 0 


= sup  < / Un{Van{x)  • (j))  + an{x)div{(f)))dx  : \(f)i{x)\  < 1 V(a;)  € D 
(j>ec^{n)^  I in 

We  choose  p,  q with  1 < p < 2 such  that  - + - = 1.  Then 

P Q 

Un  ^ u in  L^{Q)  and  ^ ct  in 

Thus  for  any  (j)  G with  |0j(a;)|  < 1 V(x)  G D we  have 

/ u{V a ■ 4>  + adiv {(/))) dx  = liminf  / ttn(Vo!„  • 0 + o;„div(0))dx 
in  in 

= liminf  / u„div(o;„</>)dx 

n->oo 

< liminf  / 0!„|Vu„|. 

rn.oo 

Hence  taking  the  supremum  over  all  (f)  of  the  left  hand  side  of  the  above  we  have 

/ q;|Vu|  < 

in 

Finally  using  this  result, 

infF(u„)  = liminf  ( / o;„|Vu„|  + {un-Uo)^dx 

Vin  in  y 

> liminf  / a„| Vu„|  + lim inf  / {un  — Uofdx 

n^oo  n-)-oo 

> / q;|Vu|  + (u  — Uo)^dx. 

in  in 


liminf  / o;„|Vu^ 

n-Fcx)  In 


The  last  inequality  follows  by  convexity  of  the  second  integral  and  by  the  fact  that 
we  can  also  choose  {«„}  to  converge  weakly  to  u in  L^(D).  Thus  F has  a minimizer 
in  L2(Q)  nBC(D).  □ 
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2.8  Other  Selective  Smoothing  Models 

We  mention  here  some  other  proposed  selective  smoothing  models  used  for 
image  restoration. 

One  such  model  was  proposed  by  Chambolle  and  Lions  [8]  which  uses  a com- 
bination of  TV  and  isotropic  diffusion.  Here  edge  preservation  is  achieved  using 
total  variation  where  the  image  gradient  is  above  a certain  threshold  and  smoothed 
isotropically  where  the  gradient  is  below  that  threshold.  The  model  is  the  minimiza- 
tion problem 

min  ^ f iVuPdx-f  f (iVnl  — |) -b  ^ f (u  — lYdx. 

BV{n)nL^{u)  2e  J\yu\<(  J\vu\>e  2 ^ Jn 

with  given  corrupted  image  I € BV  (fl)  nL^(n).  This  model  will  be  discussed  further 
in  chapter  3. 

A model  proposed  by  Chen,  Levine,  and  Stanich  [10]  is 

min  f (bix.Du) + —{u  — iV 
BV(n)nL2(Q)  2 

where 

(j){x,z)  = < 

where  ,5  > 0 and  1 < a < p{x)  < 2.  For  example  we  can  use 

^ ^ l + k\VG^*I{x)\' 

This  model  uses  intermediate  values  of  the  exponent  of  of  |Vn|  for  anisotropic  diffu- 
sion, that  is  1 < p < 2 as  well  as  TV-based  diffusion  (p  = 1)  along  edges  and  isotropic 
diffusion  in  homogeneous  regions  (p  = 2).  The  anisotropic  diffusion  is  used  for  where 
the  difference  between  noise  and  edges  is  unclear. 

Another  model  proposed  is  by  Blomgren,  Chan,  Mulet,  and  Wong  [4]  which 
avoids  the  difficulty  of  choosing  the  threshold  e or  ^5  as  used  in  the  above  models. 
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The  problem  is 

min  [ 

Jn 

where  lim|v«|-).oP  = 2 and  lim|vu|->.ooP  = 1-  This  model  also  uses  a combination  of 
edge  preserving  TV-based  diffusion  as  well  as  anisotropic  and  isotropic  smoothing. 


CHAPTER  3 

SELECTIVE  SMOOTHING  AND  PARTIAL  REGULARITY 


3.1  Introduction 


In  the  first  two  sections  of  this  chapter,  we  focus  on  the  problem 


(3.1) 


where  is  the  following  convex  function  defined  on  R" 


0 C R"  is  a bounded  domain  with  Lipschitz  boundary,  and  I G L°°(D)  Cl  BV{fl) 
is  given.  Such  a functional  for  image  restoration  was  considered  in  Chambolle  and 
Lions  [8].  Here  the  restored  image  is  taken  to  be  the  minimizer  of  a combination  of 
the  total  variation  and  the  squared  L^(Q)  norm  of  the  gradient.  That  is  we  minimize 


edge  preserving  via  TV-based  diffusion  where  edges  are  more  pronounced  (|Vu|  > e). 
Without  loss  of  generality  we  take  e = 1 as  in  (3.1). 

In  order  to  define  for  u G BV (D)  we  note  that  since  the  gradient  of  u is  a 
measure  Du,  we  can  decompose  Du  into  its  absolutely  continuous  and  singular  parts 
with  respect  to  Lebesgue  measure,  that  is 


See  Evans  and  Gariepy  ([12])  for  a complete  discussion.  Then  we  define  ([15]) 


Thus  we  expect  to  have  isotropic  diffusion  where  the  image  is  uniform  (|Vu|  < e),  and 


Du  = Vu  dx  -t-  D^u. 
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with 


[ \D^u\  = [ d\D^u\  = |Il"u|(Q). 

Ja  Jn 

It  is  important  to  note  ([30]  or  [15])  that  the  functional  J can  also  be  defined 


by 


J{u)=  sup  \~  [ ( dx  : \(f){x)\  < 1 yx  E Q. 

</>ec^(n,R")  I JO  \2  / 

Using  this  definition,  we  see  that  the  functional  J is  lower  semicontinuous  with  re- 
spect to  convergence  in  The  proof  is  similar  to  that  of  the  proof  of  lower 

semicontinuity  of  the  BV  seminorm.  Since  we  have  established  lower  semicontinu- 
ity of  J,  is  now  straightforward  to  show  that  there  is  a unique  solution  to  (3.1)  by 
standard  methods. 

We  now  state  the  two  main  regularity  results. 

Theorem  3.1.1  If  u is  the  solution  to  (3.1),  then  for  any  given  0 < /r  < 1 there 
exist  positive  constants  cq  o,nd  kq  depending  only  on  n and  /j,  such  that  if 


^ [ \Du-l\<  eo 

^r\  J Brio.) 


\Br\  JBr{a) 

holds  for  some  Br{a)  CC  and  for  some  I G R”,  with 

rC  (l  -I-  |]/]|L<x>(n))  < kq  and  |/|  < 1 - 2/r, 
for  some  constant  C depending  only  on  n and  Q then, 


D^u\{Br/2{a))  = 0 and  |Vu|  < 1 — p,  on  Br/2{o) 


and  u solves 


—Au  = I — u on  Br/2{a). 
Hence  u € C^’°‘{Br/2  (a))  for  any  a < 1. 


Theorem  3.1.2  Let  u be  as  in  Theorem  (3.1.1/  7/£"({|Vu|  < 1})  > 0,  then  there 
exists  a nonempty  open  region  E on  which  u is  |Vr/|  < 1 and  u solves 


—Au  = I — u on  E. 
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In  addition  we  have  |Vu|  > 1 a.e.  onH\E. 

It  is  actually  straightforward  to  show  that  Theorem  3.1.2  is  a direct  conse- 
quence of  Theorem  3.1.1  using  standard  results  from  measure  theory  and  analysis. 
Thus  from  Theorem  3.1.2,  we  do  indeed  have  smoothing  where  |Vu|  < 1. 

Here  we  should  point  out  that  regularity  results  were  obtained  in  Anzellotti 
and  Giaquinta  [1]  for  minimizers  in  BV (fl)  of  functionals  of  the  form 

f {F{x,  Du)  -f  G{x,  u)) 

Jn 

where  F{x,p)  is  a convex  function  in  p with  Ci|p|  < F{x,p)  < C2(l  + |p|)  for  all 
p e R"  and  G{x,  z)  satisfies  certain  continuity  conditions  in  both  x and  z.  In  our 
case,  G{x,z)  = \j2[z  — I{x))‘^  with  only  the  stated  assumption  on  I. 

The  proof  of  the  above  theorems  follows  the  ideas  of  Tonegawa  [26],  where  the 
above  theorems  were  proved  for  the  minimizer  to  the  plasticity  functional 

min  < / (fiDu)  — / fudx  > 

u€BV{n)nL^Q)  j 

for  p as  above  and  / E L°° . 

3.2  Proof  of  Theorem  3.1.1  and  Theorem  3.1.2 

First  we  will  show  that  the  solution  u to  (3.1)  is  in  L°°{D).  To  prove  this  we 
could  consider  the  time  evolution  problem  corresponding  to  (3.1),  as  in  the  case  of 
the  previous  selective  smoothing  problem  with  parameter  function  a,  prove  an  L°° 
bound  for  the  time  dependent  solution  u(x,  t),  and  then  consider  the  time  asymptotic 
limit  u.  Thus  we  would  conclude  as  in  that  case  that  u E L°°(Q).  The  next  lemma 
however  provides  a proof  of  this  without  having  to  consider  the  time  evolution  of 
(3.1). 

Lemma  3.2.1  If  u is  the  solution  to  then  u E L°°{D).  In  fact,  we  have 

||«||L°°(n)  < ll^||L°°(n)- 
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Proof:  Let  be  defined  on  R”  by 


^e{p) 


IpI  - if  bl  > 1, 


and  consider  the  minimization  problem 


min  I [ (Pf(Vu)  + - [ (u  — I)^dx\. 
uew^^{u)nLHQ)  [J^  2 Ja  j 

By  standard  methods,  there  is  a unique  solution  to  the  this  problem.  We  follow 
a standard  truncation  argument  where  we  fix  e and  f > 0 and  let  v = 

Noting  that  v G fl  with 

{Vue  if  Ue  < ^ 

0 if  > t, 

we  have 


[ ^ f {u^-lfdx<  [ (pe(Vu)  + i [ {v-lfdx, 

Jn  ^ Jn  Jn  ^ Jq 

and  thus  after  subtracting 

/ PeC^u^)dx+  / {u^  — I)‘^dx<  / {t  — I^dx. 

J {Ue>t}  J {Uc>t}  J {Ue>t} 


(3.2) 


Hence 


f (Ue  — I)^dx  < [ {t  — ly 
J {uc>t)  J {ue>i\ 


dx. 


Hue>t}  J {Ue>t) 

But  setting  t = ||/||L°°(n)  we  see  that  if  ess  sup  u^>  t then 


^dx 


which  contradicts  the  above,  hence  ess  sup  < ||/||L°°(n)-  Applying  a similar  ar- 
gument to  u = max(Uf,  — t)  for  t = ||/||L°°{n)  we  get  ess  inf  > — ||/||Loo(n)  and 
thus  |b6||L°°(n)  < ||L|U°°{n)-  Furthermore,  letting  u = 0 in  (3.2)  we  see  that  is 
bounded  in  n L'^{Q)  C BV{VL)  fl  L?(Vt)  independent  of  e.  Thus  there  is  a 
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u G BV{Q)  n and  a subsequence  of  {ue},  still  denoted  by  {ue},  such  that 

^ u strongly  in  L^(Q),  ^ u weakly  in  and  ^ u a.e  in  Q.  Letting 

e — > 0 in  (3.2),  noting  that  (/?£-></?  uniformly,  lower  semicontinuity  of  the  functional 
/ (fiyv)  defined  on  and  weak  lower  semicontinuity  of  the  second  term  on 

Jn 

the  left  hand  side,  we  get 

f yj(Vit)  + ^ f {u  — lYdx<  f ¥^(Vu)  + ^ f {v  — I^dx 
Jn  ^ Jn  Jn  2 

for  all  V G VL^’^(r2)nL^(r2).  We  now  note  ([15])  that  for  any  v G BV{Q)r\L‘^{^)  there 

exists  a sequence  in  C°°(Q)  such  that  / (p{Vvn)dx  — >•  / v?(Vu)  and  ^ u in 

Jn  Jn 

L^(Q),  and  since  v G L‘^{Q)  from  the  construction  of  ([15])  we  can  also  take  > n 
in  Z/^(fl).  Therefore  we  see  that  the  above  holds  for  all  v G BV (Q)  fl  L^(fi)  as  well. 
Hence  u solves  (3.1).  By  uniqueness,  u = u.  By  the  uniform  L°°  bound  for  and 
the  convergence  of  Uj  to  u a.e.  in  Q we  have  u G L°°{Q)  with  ||n||Loo(n)  < ][/[|L°°(n)-  D 
We  note  here  that  in  the  above  proof,  we  could  have  chosen  instead,  for 
instance,  a regularization  of  (3.1)  using  an  appropriate  smoothing  of  I with 


||-fe||L«>(n)  < PI|L«>(n) 

and  If  I in  T^(fl),  instead  of  We  also  mention  that  in  passing  to  the  regu- 
larized problem,  we  avoid  having  to  consider  the  singular  measures  in  the  truncation 
argument. 

Throughout  the  rest  of  this  section,  we  fix  ^ > 0 and  let  ci,  C2, . . . represent 
constants  that  depend  only  on  n,  /x,  u,  f2,  (/?  and  possibly  I.  Unnumbered  constants 
will  be  clearly  labeled  on  what  they  depend. 

From  Tonegawa  [26]  we  have 

Lemma  3.2.2  Let  u G BV{Br{a))  and  h G C^{Br{a))  with 

sup  |V/i|  < 1 — /X, 

Br{a) 
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then 


/ ip{Du)  — / (p{Vh)dx  >11  \D^u\  + V{u  — h)  ■ Vhdx 

J Br{a)  J Br{a)  J Br{a)  J Br{a) 

+ [ D^u-Vh  + ^ [ \Vu\dx  + l [ \\/{u-h)\^ 

J Br{a)  2 yBr(a)n{|Vu|>l}  2 ys^(a)n{| Vu|<l} 


We  fix  B2r{a)  CC  Q.  Let  u be  a Lipschitz  function  defined  on  B2r{a)  and 
assume  there  exists  an  / € R"  with  |^|  < 1 — 2/j,  such  that  sup^^^(q)  | Vu  — /|  < for 
(5  > 0 and  0 < /3  < 1 to  be  chosen  later.  Also  let  v be  defined  by  v{x)  — v[x)  — I • x. 
Let  Tje  be  the  usual  mollifier  on  R"  and  denote  = rjrp*v  and  = r]r^*v.  We  also 
have  the  estimates  ([24]) 


sup  \Vvp  — ^1  = sup  iVn^l  < , 

Br(a)  Br{a) 

sup  \vp  — v\  = sup  \vp  — v\<rP  sup  |Vu^|  < 

Br{a)  Br{a)  Br{a) 


(3.3) 


sup  |x  - yl  ^|Vu/)(x)  - Vu^{y)| 

Br{a) 

< CiT^  sup  |Vn  — /|  sup  \x'  — y'\~^\’qi{{rP)~^x')  — r]i{{rP)~^y')\ 

Sr  (a)  x'i^y' 

< C2^25^-<5  = 

Now  for  any  f G [|,r]  there  exists  a unique  solution  ([14])  w G H^{Bf[a))  fl 
C^’^{Bf{a))  with  6 G (0, 1)  for  the  problem 


—Arc  = I — w on  Bf{a) 
w = vp  on  dBf{a). 

Lemma  3.2.3  The  solution  w to  (3 A)  satisfies 

||s;||L°°(Sf(a))  < \\vy\\L°°{dBf{a))  + (^(s.,  fi) ||/||L°°(n) • 
Proof:  Let  A be  the  standard  smoothing  of  I such  that 

||A||L«>(n)  < ||^||L“°(n), 


(3.4) 


le^I 
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a.e.  in  Bf{a),  and  4 — > / in  L'^{Bf{a)).  There  exists  a smooth  solution  ([14])  to 

—Aw^  = — We  on  Bf{a)  (3.5) 

We  = vp  on  dBf{a) 

with  the  estimate  ||n;e||/,oo(B,-(a))  < ||^’^|U°°(aiJj={a))  + C'(?^,  f^)||^IU~(o)-  Let  w be  the 
solution  to  (3.4).  Since  we  can  also  bound  We  in  H^{Bf{a))  independent  of  e,  there 
exists  a subsequence  of  {we},  still  denoted  by  {we},  such  that  We  ^ w a.e.  in  Bf{a). 
Let  e — > 0 in  the  above  estimate  for  We  and  the  lemma  is  proved.  □ 

We  now  have,  using  the  above  lemma,  the  following  estimates  ([1,  14]): 


||'^||L°°(B;(a))  < \\vp\\L°°{dBf(a))  + C'(?^)  ^^)  ||Ll|L°°(n) , 

sup  r^\x  — y\~^\Vw{x)  — Vu;(y)|  + sup  |Vru  — /|  (3.6) 

Bf[a)  Bf{a) 

< + r\\I  - wWi^i^Bfia))), 


|Vtu(3:)  — Vw{y)\ 
x,yeBf/2{a)  \x-y\^ 


Note  that  from  (3.3)  we  also  have  \\I -w\\Loo(^Bf{a))  < ll^^||L°°(aSj=(a))+C'(?^,  L!)||/||Loo(n). 
Lemma  3.2.4  Suppose  there  is  a v E C^’^{B2r{a))  and  I G R”  with  |/)  < 1 — 2/u, 


sup  I Vu  — /|  < 

B2t{o) 

and  sup^2^(„)  |n|  < C„  where  Cu  is  a constant  depending  only  on  ||u||£,oo(q).  Let  vg,f, 
and  w he  as  in  the  previous  discussion.  Then  there  exists  constants  C5  and  Ce  such 
that  if  P < C5  and  r{Cu  + C{n,  fl)||/||i<x)(n))  < Cg  then  we  have 


f (p{Du)  — f ip(Vw)dx>  f {u  — d'W'  ^ 

JBfia)  JBf{a)  J dBf{a) 

f {u  — w){I  — w)dx -\- p,  f \D%\  + ^ f \Wu\dx 

J Bf{a)  JBf{a)  2 J Bf{a)n{\Vu>l} 
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I [ \V{u-w)\‘^dx 


> 


[ {u-vp)^dH^  ^ f {w-  Ifdx  [ {u-I) 

JdBf{a)  C'n  I JBf{a)  ^ J Bf  {a) 


dx 


r(“) 

+n  [ \D^u\  + f^  [ \Vu\dx+^  f 

JBf(a)  ^ ./Br(a)n{|Vu|>l}  ^ J Bi 


|V(u  — w)\'^dx 


'Bf{a)  ^ JBr(a)n{|Vu|>l}  ^ JSf(a)n{|Vu|<l} 

Proof:  From  (3.6)  and  by  assumption  on  I we  see  that 

sup  \Vw\  < sup  \Vw  - /|  + |/|  < + r\\I  - ty||L°°(B-(a)))  + 1 - 2/r 


B 


f{a) 


Bf{a) 


< + ^(lklU°°{aSf(a))  + C{n,  r2)||/||L°°(n)))  + 1 — 2/r 

^ + fiCu  + C{n,  f2)||/||i,cx>(n)))  + 1 — 2^. 

Choose  C5  and  ce  such  that  < C5  and  r{Cu  + C{n,  i7)||/||x,oo(f^))  < Ce  imply  C3(/5‘*  + 
/■(Cu + C(n,  f2)||/||L°°(n)))  < fJ>-  Thus 

sup  |Vu;|  < 1 — /i. 

5f(a) 

The  conditions  of  Lemma  (3.2.2)  now  hold  for  v = w.  Substituting  in  w for  v in  the 
inequality  in  Lemma  (3.2.2),  integrating  by  parts,  and  using  Young’s  inequality  for 
{u  — w){I  — w)  — —{u  — I){w  — /)  + (/  — the  Lemma  is  proved.  □ 

We  now  introduce  the  following  definition,  notation,  and  results  ([15]  and 

[16]). 

Definition  3.2.5  A function  u G BV{Q)  is  a local  solution  in  Q if 

[ if{Du)  + l-  [ {u-I)^dx<  [ (p{D{u  + C))  + ^ [ {{u  + Q-lfdx 
Jn  ^ Jn  Jci  ^ Jn 

for  any  f G BV{Vt).  Thus  the  solution  u to  (3.1)  is  a local  solution.  The  following 

first  variational  formula  is  from  Hardt  and  Kinderlehrer  [15]:  if  u is  a local  solution 

then 

j a -V(^dx  + I a-^|D*u|  = — j {u  — I)Qdx 

fc/  57  */  57  «/  57 


(3.7) 
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where  C is  any  function  in  BV{Cl)  with  <<  \D^u\,  ^ is  the  Radon-Nikodym 
derivative  of  with  respect  to  \D^u\,  and  a G L^{fl)  is  the  stress  tensor  defined 
by 

{(fp{Vu)  in  Qq 
D^u/\D^u\  in  fis- 

Here  D^u/\D^u\  denotes  the  Radon-Nikodym  derivative  of  D^u  with  respect  to  \D^u\ 
and  = Qa  U r^s  is  the  decomposition  of  Q with  respect  to  the  mutually  singular 
measures  £”  and  \D^u\.  Clearly  lo'(n)|  < 1.  Note  that  cx{u)  depends  only  on  the 
local  solution  u.  In  the  sequel  we  will  write  a instead  of  cr{u)  and  write  the  left  hand 
side  of  (3.7)  as 

a ■ DC,. 


/ 

Jn 


Also  note  if 


/ 

Jo. 


a ■ DC,  = — j {u  — I)C,dx 

Jn 

holds  for  arbitrary  ^ £ BV (17)  for  some  u where  a is  defined  as  above,  then  u solves 


(3.1).  In  fact,  for  arbitrary  v G BV{Q)  we  take  ( = v - u,  noting  that  by  convexity 
of  (f  we  have  (p{Vv)  — t^(Vu)  >V{v  — u)  ■ ipp{Vu)  on  f7a,  and  that  on  Dg  we  have 


" , , J 

/ \D^v\- 

Jn,  Ji 

\D^u\  > / 

rts  ^ • 

D‘‘lv  — u)  ■ — - 

i>.  \D-n\ 

Next,  define  the  functional 
= [ (p{Dv)  - f 

J Br(a)  J dBr(a 


o ■ :v  dn 

(a)  - a\ 


“ -ft 

^ J Br  (a) 


dx. 


r(a)  J dBr{a)  H ^1  ^ J Br  (a) 

We  have  a second  variational  formula  ([16])  involving  Jr,a- 

Theorem  3.2.6  The  function  u G BV{C1)  is  a local  solution  if  and  only  if 


for  every  v G BV{fl)  and  r < dist{a,dQ). 


49 


We  use  the  above  theorem,  Lemma  3.2.4,  and  estimates  (3.3)  to  obtain  the  following 
inequality  for  the  solution  u to  (3.1): 

Lemma  3.2.7  Let  v,  I be  as  in  Lemma  3.2.4  with  r(C„  + C(n,  f2)l|/||ioo(n))  < cq,  w 
as  in  (3A),  and  u a solution  to  (2>,\).  Then 

f \D^u\+  f \Vu\dx+  j \V{u  - w)\‘^dx 

< ct  f \u  — v\ d'H”'~^  + c^r"' . 

JdBf{a) 

Proof:  Clearly  a solution  u to  (3.1)  is  also  a local  solution.  By  the  second  variational 
formula  and  Lemma  (3.2.4)  we  have 

[ 0^  ■ 7“ — —Au  — VR)dTr~^=  [ cr  • -r ^Au  — w)  dW^~^ 

JdBf{a)  P ~ JdBf[a)  P “ 

> f ip{Du)  + 1:  [ {u  — lYdx  — I cp{Vw)dx  — ^ / {w  — I)^dx 

J Bf{a)  2 J Bf{a)  J Bf(a)  2 J Bf{a) 

>f  {u  - V0)^  d'H'^-^  + [ |L>*u|  + ^/  \Vu\dx 

JdBf{a)  on  JBf(a)  2 Vu|>l} 

+l  [ \V{u-v)\^dx. 

The  lemma  is  thus  proved  by  then  using  the  estimate  for  \Vw\  obtained  in  Lemma 
(3.2.4),  the  estimate  for  \v  — v^\  from  (3.3),  and  by  noting  that  |cr|  < 1.  □ 

The  proof  of  the  following  lemma  is  based  on  Hardt  and  Kinderlehrer  [16], 
and  Tonegawa  [26]  with  some  necessary  modifications. 

Lemma  3.2.8  . Suppose  u is  a local  solution  in  Q,  of  our  minimization  problem, 
B2r{a)  CC  Ll,  r{Cu  + C{n,Ll)\\I\\L°°{n))  < and  v € C^’^{B2r{a))  with 


sup  |Vu|  < 1 — pt, 

B2r{a) 


suPB2.(a)  \v\<Cu,  and 


for  all  r < p < 2r. 
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Then  there  exists  positive  constants  Cg  and  Cio  such  that  if 


^v}r\  B2r{a))  < cgr"' 


then 


\u  — v\ 


L°°(Br(a)) 


< Cio  (>C"({u  ^v}n  52r(a)))" 


Proof:  First  we  note  that  the  function  p satisfies  |p|  — A < ip{p)  < \p\  for  all  p G R", 
some  A > 0.  By  convexity  of  p we  have  cp{p)  < pp{p)  - p + <y?(0)  for  all  p G R".  Hence 
we  have 

\Du\  = |Vu|c?2;  + \D^u\  < ip{Vu)dx  + \D^u\  + \dx 
< ipp{Vu)  ■ Vudx  + \D^u\  + (A  + (f{0))dx  = a • Du  + Xdx. 


Let  0 : R -4  R be  a bounded,  increasing,  piecewise  differentiable  function 
with  9'{t)  < 1 for  almost  all  t.  Let  0 < p < h and 

1 in  Bp{a) 

r]{x)  = < {h  - {h  - \x  - a\)  in  Bh{a)\Bp{a) 

0 in  D\Bh{a). 

Now  apply  the  first  variational  formula  to  C = p6{u  — v)  to  get 


v)dx 


In  order  to  obtain  a lower  bound  for  pa  ■ D[6[u  — u)]  we  use  the  above  properties  of 
(p.  We  have  D[9{u  — u)]  = 9'{u  — v)D{u  — v)  and  hence  by  noting  the  bound  of  |Vu| 


[ \m~  - 

w)]|  < [ 9'{u  — v)\Du\+  [ 9'{u 

JBp{a) 

J Bp{a) 

J Bp{a) 

</  . 

'’{u  — v)p{Du)  + / 

(A  + 1)9' {u  — v) 

JBp{a) 

J Bp{a) 

- f ‘ 

''(u  — v)a  ■ Du  + / 

(A  + 1)9' {u  — v) 

JBp(a) 

J Bp{a) 
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j 9'{u  — v)a  ■ D{u  — v)  + j 9'{u  — v)a-Dv 
J Bp{a)  J Bp(a) 

+ f (A  + 1)9' {u  - v) 

JBp{a) 

< f rja  ■ D[9{u  — u)]  + f C\9'{u  — v) 

J Bh(a)  J Bh{a) 


' Bh{a)  JBh(a) 

for  some  constant  C\  depending  only  on  A.  Therefore 


f t]\D[9{u  - v)]\  < (h-p)  ^ f - t;)|cix  + C'A|suppr7^( 

JbM  j Bh{a)\Bp(a) 


U — V) 


' Bh{a) 

+ \\u 

For  0 < A:  < s we  choose  9 as 


\9{u  — v)\dx. 


Bh(a) 


0 for  t < A: 


9{t)  = '^  t — k lor  k < t < s 
s — k for  t > s. 

Now  let  ^(A:,  h)  = Bh  C\  {u  — v > k}.  Clearly  support  [r]9{u  — v)]  C ^(A:,  h).  Thus 
[ |T»[0(n-t;)]|<((h-p)-i  + 2||/|Uoo(n))  [ \9{u  - v)\dx  + C,\A(k,h)\ 

J Bp{a)  J Bh{a) 

By  assumption,  |t4(0,p)|  < \\Bp{a)\  for  r < p < 2r.  Thus  we  see  that 

C^{{9{u-v)  = fl}nBp{a)}  ^ 1 
\Bp[a)\  2 

Then  apply  the  isoperimetric  inequality  for  s > A:  > 0 to  get 

n—  1 

{s  - k)\A{s,p)\'^  < I [ \9{u-v)\^dx]  < Cn  [ \D[9{u-v)]\ 

\JBp(a)  J JBp{a) 

< Ci2((/i  - + ||/|U-(n))  [ \9{u  - v)\dx  + ci3\A{k,h)\. 

JBh(a) 

So  since  h < 2r  we  have 


(s  — A;)|^(s,  p)|  " < Cu(h  — p)  ^ I \9{u  — v)\dx  + Cu\A{k,h)\ 

JBh(a) 
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And  since 

f \6{u  — v)\dx  < {s  — k)\A{k,  h)\, 

we  arrive  at 

|A(s,  < cu{{h  - p)-^  + {s-  k)-^)\A{k,  h)\ 

for  every  r < p < h < 2r  and  s > A:  > 0.  We  now  apply  Lemma  2.1  ([16])  to  obtain 
the  upper  bound. 

The  lower  bound  for  n — n is  obtained  by  using  a similar  argument  for  0 < 
k < s < oo, 


0 


for  t > —k 


m = 


-t  — k for  — s < t < — A: 


s — k for  A < — s, 
and  A{k,  h)  = Bh  ^ {u  — v < —A:}.  The  lemma  then  follows.  □ 

Now  define  the  energy  function 

1 


iVnjd. 


+ 


\Br\  t^Br(i)n{|Vu|>l} 

[ \Vu-l\^dx+  [ 

J Br(x)nUVu\<l)  JErix)  J 


^Br(x)n{|Vu|<l}  J Bt{x) 

The  next  theorem  provides  a decay  estimate  for  <L. 

Theorem  3.2.9  If  u solves  (3 A)  with  Br{a)  CC  0,  A G R"  with  |/i|  < I — p,  then 
there  exist  positive  constants  cu,  e,  k,  Cis,  c\%,  and  cn  such  that 


<E>(4r,  Ai,  a)  < e 


and 


implies 


where 


r < K 


^{cor,  I2,  a)  < -<l>(4r,  h,a)  + ci^r 


\h  -h\<  Ci6‘J>(4r,  Ai,  a)  2 + ciyr. 
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The  above  theorem  (a  proof  for  a similar  result  will  be  given  in  the  next  section)  and 
subsequently  Theorems  3.1.1  and  3.1.2  now  follow.  See  for  instance  Tonegawa  [26] 
or  Hardt  and  Tonegawa  [17]. 

3.3  Partial  Regularity  for  the  p-Laplacian 

Instead  of  restricting  our  model  to  the  Laplacian  as  in  the  last  section,  we 
may  also  include  the  p-Laplacian.  So  we  can  consider  the  problem 

min  I f (fo{Du)-\ — [ (u  — I)^dx\ 

«esv{n)nL2(n)  2 j 

where  ipo  is  the  following  convex  function  defined  on  R" 


(po(^) 


where  1 < p < 2,  1/p  -h  1/q  = 1,  / e L°°{Q)  fl  and  f2  C R"  is  a bounded 

domain  with  Lipschitz  boundary.  Again  as  in  Anzellotti  and  Giaquinta  [1]  or  Hardt 
and  Kinderlehrer  [15]  we  may  define  the  above  functional  on  BV  as 


(3.8) 


) p 

if  |2;| 

1 , , 1 

if  jxj 

1 q 

/ (po{Du)  = / po{Vu)dx+  / \D‘ 

J ^ t/  V rz 


u\. 


For  applications  to  image  restoration  this  functional  is  a combination  of  anisotropic 
diffusion  (1  < p < 2)  and  TV  diffusion. 

Due  to  the  singular  nature  of  cpo  we  instead  consider  a regularized  version  of 
the  above  minimization  problem,  namely 


[ <fe{Du)  -t-  J [ {u-  Ifdx 
Ja  ^ Ja 


mm 

ueBV(n)nL2(n)  [ 

for  the  following  convex  function  pj,  e > 0,  defined  on  R"  by 

-(jrrp  -I-  if  jx]  < 1 

(p,{x)  = { ^ I 

(1  + e)P/^~La:| (1-l-e)^/^  ^(p— 1 — e)  if  jx]  > 1, 

P 


(3.9) 
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where  $7  C R",  I,  and  p,  are  as  before.  As  above, 


We  note  that  for  e > 0,  (/?£  G and  is  on  the  interior  of  the  unit  ball  Ri(0). 

First  we  show  that  / p^{Du)  is  lower  semicontinuous  in  L^{Q)  for  any  e > 0. 

Jn 

Lemma  3.3.1  For  any  e > 0,  the  functional  / Pe{Du)  is  lower  semicontinuous  in 

Jn 

L^{n). 

Proof:  Let  V = {(p  E C'o(fl,R")  : \(p{x)\  < 1 Vx  G f7}.  Without  loss  of  generality 
we  can  adjust  if  necessary  so  that  <Pe(0)  = 0 and  <Pe{x)  = \x\  — K where  lx]  > 1 
for  some  constant  K.  From  Eckland  and  Temam  [11]  we  have  for  each  x G R", 

<p,{x)  = sup{a:  ■ y - ip*{y)  : y G R",  \y\  < 1} 

where  p*  is  the  convex  conjugate  of  pe  defined  by 

p*{y)  = sup{a;  • y - p^{x)  : x G R"}. 

For  the  special  case  of  e = 0 we  in  fact  have 

Pq{x)  = sup{2:  • y - \y\^/q)  : y G R",  ly|  < !}• 

The  linear  growth  property  of  actually  gives  finite  values  for  p*{y)  only  when 
|y|  < 1 as  well  as  the  fact  that 

(p*(y)  = sup{x  • y - p^{x)  : x G R",  |a:|  < 1}. 

We  note  that  we  actually  have  for  e > 0,  <p*(y)  = x*{y)  ■ y — (p^{x*{y))  where  x*{y)  is 
a continuous  function  of  y.  From  Brezis  [5]  we  have  for  any  g G 


(3.10) 
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Now  define  the  following  functional  on  BV{Q): 

J(m)  = sup  < — / u div 4>  + dx> 

0€V  [ JQ  J 


= sup / Vu  ■ 4>  — ip*{(j))  dx / 
<pev  {Jn  Jn 


^■D^u 


where  the  last  equality  follows  from  integration  by  parts.  From  the  above  discussion, 
we  easily  see  that  for  every  (j)  eV 

I Vu  ■ (j)  — (p*{(j))  dx  + f (pD^u  < j (ft{Du) 

Jn  Jn  Jn 

giving  J(u)<  / (p,{Du). 

Jn 

For  the  reverse  inequality  we  follow,  for  example,  Chen  et.al.  [10],  noting  the 
continuity  of  cp*.  Fix  e > 0.  For  any  u G BV{Cl)  there  exists  an  open  set  such  that 
support(D^u)  C and  \Of\  < e.  We  can  also  find  a 0i  G Co(fl,R")  with  |(^i|  < 1 


and 


[ D^u  ■ (j)i  > f \D^u\dx  — 
Jn  Jn 


(3.11) 


from  the  definition  of  the  TV  norm.  By  (3.10)  there  exist  a (/>2  G C'd(fl,R")  with 
\(f>2\  < 1 such  that 


/ Vu  • 02  - > (p^{Vu)d: 

Jn  Jn 


X — e. 


(3.12) 


Now  define 


01  on 


0 = 


02  on  fl  — Oe, 

Let  Tja  be  the  standard  mollifier  on  R"  and  let  0^  = T]a*<p-  Note  that  0^  has  compact 
support  for  sufficiently  small  a.  Then 

J{u)>  / Vu  ■ 0a  - <p*(0a)dx  + / 0„  • 

Jn  Jn 

Letting  a — >•  0 in  the  above  inequality  we  then  have 

J{u)  > I Vu  ■ (f)  — ip*{^)dx  + f <p  ■ D^u. 

Jn  Jn 
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> f \/u  ■ (f)2  - ip*{(j)2)dx  + f (j)i  ■ D%  - fi{e) 

Jn  Jn 

> / ip^(Vu)dx+  / \D^u\  — fj,{e)  — 2e 

Jq  Jn 


where 

/^(e)  = [ \^u\dx+\\(t>*\\L°o[Br(0)}\Oe\. 

JOe 

Clearly  /r(e)  ->  0 as  e ->  0.  The  reverse  inequality  is  now  proved.  Lower  semiconti- 
nuity then  easily  follows  as  in  Theorem  2.2.6.  □ 

We  also  have  an  approximation  lemma  as  in  Theorem  2.2.7. 


Lemma  3.3.2  Let  u G BV{fl)  n L‘^{Q).  Then  for  any  e>  0,  there  exists  a sequence 
of  functions  {un}  C BV{fl)  fi  L^(f2)  n C°°{Ll)  such  that 


Uji  ^ u in  and 

/ (fe{Dun)dx^  / ipe{Du). 

Jn  Jn 


Proof:  Fix  e > 0 and  for  simplicity  write  as  (/?.  Consider  the  function  <Pn{x) 
ip{x)  + \x\/n.  Then  ([!])  there  exists  a G BV{Vt)  fl  C°°{Q)  such  that 


1.  \\Un  - u\\l2{q)  < 1/ 


n 


2. 


3. 


f y/l  + {DunYdx  - f v^l  -I-  {Du)‘ 
Jn  Jn 

/ ipn{DUn)dx-  / (Pn{Du) 

Jn  Jn 


< 1/n 


< 1/n. 


The  proposition  there  is  actually  stated  for  functions  in  BV,  but  an  adjustment  of 
the  proof  on  which  this  proposition  is  based  ([2])  as  in  Remark  2.2.8  gives  us  estimate 
1 for  n G BV{Ll)  n Then 

/ (f{Du)  - / ip{DUn)dx  < / <Pn{Du)  - / iPn{DUn)d: 

Jn  Jn  Jn  Jn 

+ l/ni  f \Dun\dx  + I \Dt 
\Jn  Jn 
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<l/n  + l/n  \Dun\dx  + L jOui)  . 


From  2,  / \Dun\dx  is  bounded.  Letting  n 0 proves  the  theorem.  □ 

Jn 

These  lemmas  then  imply 

Theorem  3.3.3  For  any  e > 0,  there  exists  a unique  solution  u G L°°{Q.)  to  problem 
3.9.  In  fact  we  have  ||u||/,<x>(n)  < ||/||L«>(n)- 

Proof:  The  proof  is  standard  using  lower  semicontinuity  and  convexity.  The  L°° 
bound  for  u follows  as  in  Lemma  3.2.1  using  both  of  the  above  lemmas.  □ 

Now  we  state  the  regularity  theorem. 

Theorem  3.3.4  If  u solves  (S.9)  for  e > 0 and  £”({|Vu|  < 1})  > 0,  then  there 
exists  a nonempty  open  region  on  which  u is  |Vrt|  < 1,  and  u solves 

— div{(pp{Vu))  = I — u on  Cl. 

In  addition  we  have  |Vu|  > 1 a.e.  on  Cl\Cl. 

For  simplicity  we  consider  the  case  where  e = 1 and  we  let  (p  = (pi.  To  prove  the 
above  theorems,  we  use  a method  similar  to  that  in  Hardt  and  Tonegawa  [17]  for 
proving  partial  regularity  for  weak  solutions  u G L‘^{[0,oo],  BV{Cl))  to  problems  of 
the  form 


du 

dt 


= div^p(Vu) 


on  C or  where  is  a convex  linear  growth  function  satisfying  local  ellipticity 
and  continuity  assumptions. 

The  essential  part  of  this  result  is  Theorem  3.3.7,  but  first  we  need  some 
preliminary  lemmas.  We  have  ([17]) 

Lemma  3.3.5  Let  u G BV{E)  with  open  region  E CC  LI  with  smooth  bound- 
ary.Then  there  exist  constants  Ci,C2  < 1/2  such  that  if  p E Ri(0)  and  h G C^{E) 
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with  sup^;  \Vh  - p\  < Cia,  then  for  any  vector  pi  € -Bcia(p); 

f p{Du)  - f (p{Vh)dx-  f (pii{pi)  -Vh-  D{u-  h)  (3.13) 

Je  Je  Je 

[ {pii{pi)  ■ Pi  - Pi{pi))  ■ D{u- h)  + sup  uj{\Vh  - Pi\)^  \Vh-pifdx 

J E ^ d E 


+ 


> C2  ( / \Du 

I Er\{Du^B„(p)} 


P\+  [ 
J El 


\D{u-h)\'^  , 


' En{DueB„{p)} 

where  ci;  : R — > R is  a non-decreasing,  nonnegative  function  with  limt^o^(0  = 0- 
Recall  that 


L 


En{DueB„(p)} 


\D{u  - h)p  -I-  [ \Du-p\ 

J En{Du^B^{p)} 


means 


f ~ h)\‘^dx  -I-  f 

JEniVuGB^m)  J E\ 


I Vu 


p1+  / 

JE 


' En{VuGBa{p)}  J En{Vu^B^{p)} 

Throughout  the  rest  of  this  discussion  u will  be  the  solution  to  problem  (3.9). 

Lemma  3.3.6  Let  E CC  Ll  be  an  open  elliptical  region.  Now  suppose  h E C^{E) 
satisfies  sup^;  | Vh  — p|  < Cicr  and 

d'^h 


PhijiPi) 


dxidxj 


I — h on  E 


(3.14) 


h = vp  on  dE 


for  some  pi  E Rci<t(p)  smooth  vp.  Then  we  have 

f \Du-p\+  [ \D{n-h)\^ 

J En{Du^B^(j)}  J En{DueBa{p)} 


< C3 


( f \u  — vp\  dTE  ^ + supo;(|V/i  — pi\Y  / |V/i  — pi\^dx\  . 
\JdE  E J E / 


Proof:  Such  a solution  h exists  by  elementary  linear  theory.  A simple  modification 
of  the  proposition  in  Hardt  and  Kinderlehrer  [16]  yields  the  variational  formula 


Je{u)  < Je{v) 
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for  the  minimizer  u and  for  all  v G BV (f2)  where  Je{v)  is  defined  for  all  v G BV (f2)  fl 
L^{^)  by 

Je{v)  = [ - [ a ■ hv  + ]r  f {v  ~ I^dx, 

Je  JdE  ^ J E 

with  stress  tensor  a defined  as  in  the  previous  section.  Using  this,  integrating  by 
parts,  and  using  Young’s  inequality  for  {u  — h){I  — h)  = —{u  — I){h  — /)  + (/  — 
the  lemma  is  proved.  □ 

The  energy  functional  for  this  problem  is  given  by 
$(r,p,  a)  = -^  ( / \Du-p\+  \Du-pi 

\^r\  \J Br{a)n{DuiB„{p)}  J Br{a)n{Du€B^{p)} 

Our  next  goal  is  to  prove 

Theorem  3.3.7  There  exist  constants  Ci  and  62  depending  on  p and  p such  that 
«/ $(ro,Pi, xo)  < Cl  for  some  Tq,  xq  E Q and  some  p\  G B(:^{p),  then  there  exists 
p G B^{p)  such  that  limr_>oo  = 0.  Furthermore,  p = Du{xq). 

From  now  on  we  drop  the  ”a”  in  Br{a)  and  in  $ unless  noted  otherwise.  In  order  to 
prove  Theorem  3.3.7  we  will  obtain  estimates  for  / \u  — vp  \ dTF~^,  \Vh  — pi\,  and 

JdE 

/ I V/i  — pipdx  for  a suitable  Lipschitz  function  v.  These  will  then  be  used  to  prove 
Je 

a decay  estimate  for  $ on  a smaller  ball  and  a different  p,  whose  difference  can  be 
estimated.  This  decay  estimate  will  be  used  to  prove  the  theorem. 

Now  choose  a Lipschitz  function  v defined  on  Brj2[o)  ([26]).  From  the  proof 
of  Proposition  3 there  we  have 

Lemma  3.3.8  If  5^/2  CC  Q,  and  for  v defined  as  above,  then  we  have 

sup  I Vu  — Pi|  < 

Bt/2 

with  S such  that  (1  - = 1 + that  is  5 = 777-^—77.  In  addition  we  have 

n 2n  8(n  + l) 

C^{Br/2  F{u^  t;})  < C4r”$(r,pi)^"‘‘‘^. 
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From  the  above  lemma  we  get  an  L°°  estimate  for  u — u. 


Lemma  3.3.9  If  \pi  — p|  < Cicr,  then  there  exists  an  e > 0 such  that  $(r,pi)  < e 
implies 


sup  I Vn  — Pi  I < $(r,  pi) 


25 


B 


r/2 


for  V defined  as  above  with  5 = 
ciently  small: 


1 


8(n  + 1) 


and  with  the  following  estimate  for  r suffi- 


/ 1 \ 1/n 

^5  ( f Br/2] 

C5^(£"({u^r;}n5,/2)'/". 


Proof:  This  follows  from  the  estimate  of  £”(5^/2  Fl  7^  n})  in  Lemma  3.3.8,  the 
bound  |Vn|  < |p|  + cicr  + <f>(r,  pi)^^,  and  from  Lemma  3.2.8  where  we  use  the  above 
bound  for  |Vw|  in  that  proof  instead  of  the  bound  |V?;|  < 1 used  there.  □ 

We  can  now  estimate  |V/i  — pi|  on  E for  any  f < r j2  for  the  solution  h to 
(3.14).  Let  be  the  smoothing  of  v as  described  in  the  previous  section  with  all  the 
accompanying  estimates.  Then  from  the  linear  theory  ([14])  as  in  the  last  section, 
taking  /?  = $(r,pi),  we  have  ||/i||l°°(£:)  < Ce  (see  Lemma  3.2.3)  and 


sup|V/i-pi|  < C7(4>(r,pi)‘^ + r||/i- /||loc(£;))  < C8($(r,pi)'^  + r). 

E 

Let  T : fl  ^ R”  be  an  appropriate  transformation  such  that  \ih  = ho 
u = uo  T“\  vg  = v/3  o T~^,  i = I o E is  an  ellipse  centered  at  a,  T{E)  = B'  is 
a ball  of  radius  f centered  at  T{a)  and 


— Ah  = I — h on  B' 
h = vg  on  dB' . 


(3.15) 


We  note  that  the  Jacobian  of  T,  which  depends  only  on  the  eigenvalues  of  the  matrix 
[Pii/;(Pi)j)  is  bounded  from  above  and  away  from  0.  Furthermore  these  bounds  can 
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be  made  independent  of  a and  pi  due  to  uniform  ellipticity.  We  have  ([14]) 
\Vh{x)-Vh{y)\  ^ C9 


sup 

B 


< 


\x 


yjl/2  ^n+1/2 


[ \v0\dn^-^  + f\\h-i\\L. 

JdB' 


where  B C B'  is  concentric  with  B'  with  radius  f/2. 

Choose  E C Brf4  centered  at  a with  diamE  = r/8  and  such  that  both 


[ \u  — v\  d'H"  ^ < — [ \u  — v\dx, 
JdE  f Jb^/4^ 


and 


[ \u- ub.^  - Pi  ■ {x  - a)\d'W  ^ < — [ \u  - ub.^  - Pi  ■ {x  - a)\dx  (3.16) 

JdE  f JB^/i 

hold,  where  Ub^^  denotes  the  average  of  u over  Sr/4- 

Lemma  3.3.10  Let  v be  as  in  Lemma  3.3.8  with  its  smoothing  vg.  Then 

[ \v^  - u\  dU^-^  < Cur”(<h(r,pi)i+25  ^ $(r,pi)i+i/(2")). 

JdE 

Proof:  By  the  properties  of  vp  we  have  sup^;  j?;^  — '^1  < r$(r,pi)^+^'^,  and  from 
Lemma  3.3.9  with  the  above  integral  boundary  for  estimate  for  \u  — v\ 


f \u  — v\dB‘^  ^ < Ciir"$(r,pi) 
JdE 


l+l/(2n) 


Hence  the  lemma  is  proved.  □ 

Lemma  3.3.11  If  h solves  (3.14J  and  \pi  —p\<  cia  then 

f \Vh  - pil'^dx  < C2o(r”+^||/i  - i'||L-(E)  + r”4>(r,pi)). 

J E 

Proof:  Multiply  (3.14)  by  h — vp,  integrate  by  parts,  use  the  fact  that  (3.14)  is  a 
linear  equation  with  constant  coefficients,  and  then  use  Young’s  inequality  to  arrive 
at 

/ |V/i  — pipdx  < Ci2  / {h  — I){h  — vp)dx  + Ci3  / \Vvp  — pi\^dx. 

Je  Je  Je 
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By  the  uniform  bounds  of  V/i  and  we  see  that 


\\h-vp\\L  °°{E)  < cu{diam  E)  < C\^r. 

For  other  part  of  the  estimate  we  use  ([24]) 

/ I Vu/3  — < Ci6  / \Vv—pi\^dx. 

Finally,  to  estimate  / |Vu  —pi^dx,  the  construction  of  v gives 

J B^/2 

/ \Vv -pi\^dx  <cn  I |-0^i-PiP  + Ci8$(r,pi)  < Ci9<h(r,pi). 

Jb^!2  J B^/2^{B>u^B^{p)} 

Combining  the  above  estimates  proves  the  lemma.  □ 

We  now  arrive  at  our  decay  estimate. 

Theorem  3.3.12  There  exist  positive  constants  e,  C20,  and  k depending  only  on  n, 
n,  and  u such  that  if  ^{r,pi)  < e and  r < C20,  then  there  exists  p2  € R"  such  that 
^{Kr,p2)  < ^^{r,pi)  + C3ir  and  \p2  - pi\  < C35$(r,pi)^/^  + C34r. 

Proof:  Using  Lemma  3.3.6  and  the  estimates  obtained  in  Lemmas  3.3.10  and  3.3.11 
we  have 

f \Du-p\+  f \D{u-h)\^ 

J En{Du^B^{p)}  J En{DueB„{p)} 

< C21  +o;(c8($(r,pi)''  + r))2(r"+^||/i-/||L-(E)  +r'^$(r,pi))  . 

Letting  p2  = Vh{a),  we  now  estimate  sup  |Vh(2;)  — Vh(a)|  over  a ball  B^r  C 
E C E.  Since  h — ub^^^  — p\  • {x  — a)  also  satisfies  PDF  (3.15)  we  see 
\Vh{x) -Vh{y)\  ^ C22 


sup 

B 


-y\m 


< 


f \v~y 

JdB' 


JdB' 

+C22r^/^||h  - /||l-. 


UB^/^-pi-{x-a)\dH"‘  (3.17) 


Recall  that  the  ball  B is  concentric  with  B'  and  with  radius  r/2  and  center  T{a). 
Changing  back  to  the  original  variables,  using  (3.16),  and  finally  Poincare’s  inequality, 
the  right  side  of  the  above  can  be  estimated  as 
C23 


< 


r.ri+1/2 


f 

JdE 


ub,,,  -Pi-  {x  -a)\ dTE  ^ + C22r^^^||h  - I\\l 


63 


^ 1^^  - ^1  - pi  • (^  - 

+C24r-^/^||/i  - I\\l°-{e) 

^:;:^[  [ IDu-pildx  + c^^r^^^h- i\\lc.^e)- 

Let  E = T~^{B)  and  restrict  k as  necessary  so  that  5«;r  C E C E.  Hence  we  easily 
see,  after  changing  variables  in  the  left  hand  side  of  (3.17),  using  Lemma  3.3.10, 

^1/2  f ^1/2  r 

sup  |Vh(a:)  — V/i(y)|  < C2e( — — / \v ^ — u\ dl-E~^ -\ — / \Du  — pi\dx 

JdE 

+ K^^^r\\h-  I\\loo(^e)) 

f'  Jb,/^ 

We  thus  obtain 

sup \Vh{x)  - Vh{y)\  < C2^K^I'^{^{r,pif^^^  + 

Bkt 

+ ^{r,piYl‘^ -^r\\h~  I\\l^^e))- 

This  is  our  desired  estimate  for  |Vh(x)  — Vh{y)\. 

Now  use  the  inequality  | Vu  — Vhp  > ^|Vu  — P2p  — |V/i  — P2p  we  arrive  at 

[ \Du-p\+  f |Du-p2p 

J BKrn{Du^B„{p)}  J B^r<^{DueB„(p)} 

< C29(r"$(r,pi)^+^/(^'“)  +ta(c8(<f>(r,pi)'^  + r))V"‘^^||/i  - /||l«>(e)) 

+C29(o;(c8($(r,pi)‘^ + r))^r"$(r,pi)  + f \Vh  - p2\‘^dx). 

J B^r 

For  the  rest  of  the  proof  we  denote  $(r, pi)  by  $.  Using  the  estimate  for 
\Vh{x)  — Vh(p)|,  recalling  that  p2  = Vh{a),  and  dividing  the  above  inequality  by 
we  have 


$(ACr,P2)  < C30K  "$^+V(2n)  "o;(c8'h'^  + C8r)^<h  + C3oK<L 

+ 030(^(08$*^  + C8r)^K“"  + Kr)r\\h  - /||l°°{e). 


64 


Restrict  k again  so  that  k < . Then  restrict  <4>  and  r so  that  + 

4c3o 

C3o«:“"'a;(c8<4''*  + c^rY  < 1/4.  This  proves  the  decay  estimate  for  $.  Finally  we  derive 
the  estimate  for  \p2  — Pi|-  From  the  linear  theory  ([14])  as  applied  to  h we  have 

\P2  -Pil  = |V/i(a)  -pij  < 0327^7  / \vp  -pi  • (x  - a)  -UB'\dW~^ 

\^\  JdE 

+C33f||h  — 

< C32t47  [ \v^  - u\drr~'^  + ^32]^  I \u-pi-  {x-a)  - UE\d'H'^~'^ 

\^\  JdE  rI  JdE 

+C34r||/l  — 

Then  using,  the  boundary  estimates,  Poincare’s  inequality,  and  Holder’s  inequality 
we  get 

|P2  -Pil  < C35$(r,pi)^/^  + C34r||/l-  7||  LOO.  □ 

By  using  this  decay  estimate  iteratively,  we  then  have  Theorem  3.3.7,  as  in 
Tonegawa  [26].  We  actually  have  Theorem  3.3.7  holding  for  all  x G 5^/2 (2^0)  if 
$(ro,pi,xo)  is  sufficiently  small  by  noting  that  d>(r/2,pi,x)  < 2"<f>(r,pi,xo)  for  all 
X e Br/2{xo).  Theorem  3.3.4  then  follows  ([17],  [26],  or  [1]). 


CHAPTER  4 

A QUESTION  FOR  FURTHER  STUDY 


In  the  previous  chapter,  we  mentioned  the  partial  regularity  problem  for  the 
(unregularized)  p-Laplacian.  For  convenience  we  restate  the  problem  here.  We  ex- 
amine the  solution  to 


min  I [ ipo{Du)  + - [ (u  — I)^dx\ 

ueBV{n)nL^n)  [J^  2 ' j 

with  <po  being  the  convex  function  defined  on  R”  by 


f -|a;|P 

f P 

if  |x| 

1 , , 1 

if  |a;| 

— 

1 q 

(4.1) 


where  1 < p < 2,  and  1/p  -h  1/q  = 1.  We  wish  to  obtain  a partial  regularity  result 
similar  to  Theorems  3.1.2  and  3.3.4.  We  note  that  the  method  used  for  the  regularized 
p-Laplacian  would  not  give  regularity  on  the  set  {|Vu|  < 1}  for  the  minimizer  u,  but 
only  on  {0  < |Vu|  < 1},  as  po  is  singular  where  Vu  = 0.  An  alternative  approach 
would  be  to  apply  the  method  used  for  the  Laplacian  in  the  first  section  of  chapter 
3.  A lemma  similar  to  Lemma  3.2.2  can  be  obtained.  In  fact  we  have 


Theorem  4.0.13  Let  u 6 BV{Bj.{a))  and  h G C^{Br{a))  with  supg^(^)  |V/i|^  ^ < 
1 — p,.  Then 


f Po{Du)  - f 

J Brio)  j Br[a) 


Po{Vh)dx  > p 


j 

JBr{a 

I 

J Br(a 


\D^u\  + 


[ V(u  - h) 

J Br{a) 


V/i 

|V/l|2-P 


dx 


Vh 


+fih)  [ 

J Br(a 

/ 

J Br(a 


{Vuldx 


(a)n{|Vu|>i} 

I V(u  — h)\‘^dx. 


{a)n{|V«|<l} 
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for  some  increasing  function  f defined  on  [0, 1]  such  that  /(//)>  0 for  ^ 0 and  a 
constant  c depending  only  on  p. 


Proof:  By  definition,  where  |Vii|  > 1 we  have 


iPq{Vu)  — ipo{Vh)  — V{u  — h) 


Vh 


|Vh|2-p 


|Vu|-^ — ^--|Vh|P- V(M-h) 


p p 


|Vh|2-p 


iv^i  - ^ + ivhr 


p p 


|Vh|2-p 


> |Vm|  - + ?~~\Vh\P  - |Vu||Vhr\ 

j)  — 2)  1 

Now  let  g{x)  = a 1 — — ax^~^  for  a = |Vn|  and  x G [0, 1].  Then  since 


P 


P 


a > 1, 


g{x)  — a{l  — x^  ^)  — a{- — — ^x^) 
> a{-  + ^^xP  -xP-^). 


P 


1 J)  — 1 

Letting  g{x)  = a{-  H x^  — x^^^)  we  easily  see  that  g is  strictly  decreasing  on 

[0, 1]  with  ^(0)  = 1/p  and  ^(1)  = 0.  Now  let  /(/r)  = ^(1  — /r).  Where  |Vu|  < 1 there 
holds 


iPi{Vu)  — <^c(Vh)  — V(u  — h) 


Vh 


> c|V(m  — h)|' 


(y/|Vh|2  + e)2-P 

where  is  defined  as  in  (3.9)  and  c depends  only  on  p.  This  follows  from  the 
properties  of  the  restriction  |Vu|  < 1,  and  the  assumption  of  Vh.  Letting 

e — > oo  we  get 


yjo(Vn)  - cpoi'^h)  - V(u  - h) 


V/i 


|Vh|2-P 


> c|V(n  — /i)p 


For  the  singular  measure  we  see 


/ |D-«|  > / D‘u--^+f  |D*«|(l-|Vftr') 

J Br{a)  J Br[a)  J Br{a) 


> 


f 

J Br  [d 


(a)  |V/l|2-P 


p f \D^u\. 

J Br{a) 


Combining  the  above  results  proves  the  theorem.  □ 
However,  in  considering  the  auxiliary  PDE 
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-div( 


Vw 


|Vu;| 


2-p' 


W 


I — w on  Bf{a) 
V0  on  dBf{a). 


(4.2) 


sufficient  estimates  such  as  (3.6)  are  difficult  to  obtain.  There  are  regularity  results 
([28,  27,  19])  for  the  p Laplacian  with  1 < p < 2,  but  it  is  not  clear  if  the  bounds 
provided  are  sufficient  for  obtaining  the  decay  estimate. 
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